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SUMMARY 


This document describes performance requirements and design characteris- 
tics of a microwave scatterometer wind sensor for measuring surface winds over 
the oceans on a global basis. A review of microwave scatterometer aircraft 
and satellite measurements on the ocean is included, showing the sensitivity 
of this instrument as a remote wind sensor. This wind sensor will be used on 
the SeaSat A satellite along with several other remote sensors: a radar altim- 

eter, an imaging radar, and an optical radios ater. Together, this complement 
of instruments will provide the SeaSat A User community with scientific and 
operational type data in the areas of oceanography , meteorology and geodesy. 

The scatterometer specifications are developed from User requirements of 
wind vector measurement range and accuracy, swath width, resolution cell size 
and measurement grid spacing. The manner of transforming wind speed accuracy 
into instrument accuracy requirements is described in detai.i . Since the 
scatterometer is sensitive to both wind direction and wind speed, special tech- 
niques are developed to determine both parameters. 

From several options available, the fan beam scatterometer design has 
been chosen by the SeaSat study teams as the one most suitable when considering 
attitude control, spacecraft interface problems, and performance characteris- 
tics. A detailed analysis is performed for a baseline scatterometer design. 

The analysis includes derivation of the instrument design and its performance 
capabilities for meeting the SeaSat A User requirements. Various modes of 
operation are discussed which will allow the resolution of questions concerning 
the effects of sea state on the scatterometer wind sensing ability and to 
verify design boundaries of the instrument. 

INTRODUCTION 

Global measurement of ocean winds has been established by the SeaSat A 
User panel as one requirement of the oceanographic satellite program. Surface 
winds are required as inputs to wave forecast models and can also be helpful 


in weather forecasting. The lack of sufficient wind and pressure data over the 
oceans has thus far precluded better long range weather forecasting for 
continental areas. 

The remote sensor selected for this surface wind measurement is the 
microwave scatterometer. Aircraft and satellite scatterometer data show that, 
for windspeeds less than 20 m/sec and incidence angles beyond 25°, the ocean 
radar scattering coefficient (a°) changes about one-half dB per 10 log m/sec 
change in wind speed. 'While some authors argue whether the rate is l/2 or 1/4 
dB/10 log m/sec, the sensitivity of 0° to wind speed is still sufficient 
(even to gale-force winds) to make the scatterometer a viable remote sensor of 
wind speed and direction if proper attention is given to instrument accuracy. 

During the SeaSat . Phase A feasibility studies, several instrument designs 
were considered so that the payload could be optimally configured. The total 
instrument complement included not only a scatterometer but a radar altimeter, 
an imaging radar, and both optical and microwave radiometers. 

As scatterometer sponsor, Langley submitted to the study teams several 
candidate designs for the scatterometer which would satisfy the requirements 
set forth by the User panel (ref. l). During the August 19 73~ August 1974 
time period the Phase B Study Team, User Working Group and Instrument Working 
Group considered these options, their inpact on the spacecraft configuration, 
and their effect on User requirements. From this effort came the final 
recommended payload design which included the fan beam concept for the 
scatterometer. 

The scatterometer characteristics developed in this report are based on 
User requirements established by the SeaSat User Working Group (ref. 1). In 
an effort to obtain further insight, a specific scatterometer design with 
supporting analysis is presented. Certain design options using the fan beam 
concept are pointed out such as a high resolution mode for areas where strong 
wind speed gradients exist. 
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SYMBOLS 


A arbitrary constant 

2 

A q ground area of doppler cell, km 

doppler bandwidth across doppler cell, Hz 

O 

C speed of light, 2.996 x 10° m/s 

C element excitation distribution 

n 

d horn spacing 

D variable controlled by wind speed 

E($) electric field intensity 

fjj 'doppler shift at center of doppler cell, Hz 

f^ transmitter frequency, GHz 

F receiver preamplifier noise figure 

F($) element pattern 

G gain of antenna to center of doppler cell 

G q peak antenna gain 

H(f) transfer function of a gated mean integrator 

k wave number of ocean 

K Boltzman's constant 1.38 x 10 -2 ^ j/°K 

Kp ratio of standard deviation (o^) of 0° measurement to true value of 

a° 

IC normalized standard deviation after averaging measurements from 

pci 

multiple data cells 

Kj. ratio of noise integration period to signal plus noise integration 

period 

1 Q distance on Earth's surface from center to center of adjacent 

doppler cells along the antenna beam, km 
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L 
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** 
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N 
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R. 
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c 


length of doppler cell "base, km 
miscellaneous system loss (2 way) 
receiver loss 
integer 0, +i, +2 
integer 

receiver noise power , watts 

number of antenna elements 

numbe* of antennas used to make measurements 

number of polarizations used per antenna 

power spectral density of a random signal in bandwidth B c (one sided) 

power density of the transmitted signal at the surface of the earth 

estimate of the mean noise power 

true value of return signal power, watts 

estimate of the mean signal plus noise power 

estimate of return signal power, watts 

transmitter power, watts 

equivalent resolution size, km 

radius of the Earth, km 

distance from Earth center to satellite , km 

slant range from satellite to center of doppler cell, km 

ground range from sub~satellit e point to center of doppler cell, km 

slot spacing 

received signal power, watts 

ground spacing between center of doppler cell of front antenna beam 
and the corresponding doppler cell center of the aft antenna beam, 
km 
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surface distance separating the two orthogonal data cell measurements 
of the same incidence angle 

spectrum at the low pass filter (LPF) output 

spacing along satellite subtract between successive scans of the same 
antenna 

wave number spectrum 

time gate during each repetition period for data readout, sec 

measurement period per antenna per polarization, sec 

pulse repetition period of transmit pulses, sec 

antenna brightness temperature, °K 

reference temperature, °K 

system noise temperature, °K 

earth's surface velocity at any latitude 

earth's surface velocity at the equator 

satellite orbit velocity, km/s 

satellite ground velocity, tou/s 

component of earth's surface velocity r. the plane of the antenna 
beam 

vector sum of satellite velocity component and earth's velocity 
component in the plane of the antenna beam 

wind speed, m/s 
coordinates 

inclination of the satellite orbit 
3 cLB beamwidth of fan beam antenna in wide dimension 
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earth center angle measured from subsatellite to center of doppler 
cell, degrees 

angle between constant doppler line and line normal to satellite 
track (doppler line rotation angle), degrees 

amount of spacecraft or antenna rotation to compensate for earth’s 
velocity 

phase between adjacent antenna elements 
antenna efficiency 

nadir angle from satellite to center of doppler cell, degrees 

incidence angle of center of doppler cell, degrees 

nadir angle from satellite to direction of peak antenna gain, 
degrees 

wavelength of transmitter frequency 
ocean wavelength 
waveguide wavelength 
wind speed power coefficient 
normalised radar scattering coefficient 
estimate of a 0 

standard deviation of 0° measurement 

standard deviation of noise measurement 

standard deviation of signal plus noise measurement 

integration period 

transmit pulse length, sec, 

integration period for signal plus noise,; sec, 
integration time for noise only, sec, 

3 rlB fan beam antenna beamwidth (narrow dimension) 


<j> L latitude of the satellite 

4> rri latitude of doppler cell n 

A(f> T ^ latitude difference between the satellite and a doppler cell 

d) angular direction front antenna data cell measurement moves at the 

m 

time orthogonal mf asurement by the aft antenna is made - angle 
measured from satellite subtrack 

cj) 1 surface angle from satellite track to center of antenna beam, degrees 

$ polar angular coordinate 

angle between the earth's velocity and the plane of the antenna beam - 
ascending portion of the orbit 

ijjp angle between the earth’s velocity and the plane of the antenna beam - 

descending portion of the orbit 

Oh angle between the earth’s surface velocity and the satellite negative 

velocity - measured in the plane of the antennas parallel to the 
surface of the earth 

longitude of the satellite 

longitude of doppler cell n 

longitude difference between the satellite and a doppler cell 
Subscripts 


^Ln 

4S2 l : 


jn 

General 


HH 1st letter is transmit polarization and 

HV 2nd letter is receive polarization 

1 lower edge of doppler cell 

m integer 

n integer 

u upper edge of doppler cell 

W 1st letter is transmit polarization and 

VH '2nd letter is receive polarization 
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BACKGROUND 


Backscatter Theories - A Prospectus 

The theory of radar backscatter from rough surfaces, and particularly 

2 

from the sea, has been the subject of much recent research. Most modem 

theories describe the ocean surface by a model that involves two scales , one 

for that portion of the ocean spectrum having wavelengths much larger than the 

radar wavelength (gravity waves) and the other for the part of the ocean spec- 

trum comparable with the radar wavelength {capillary waves). Wright and 
it 

Bass and Bocharov used the small perturbation method to describe the baek- 

5 

scatter; later, Wright considered the effect of the perturbed surface lying 

on tilted planes caused by the larger structures. A similar approach was used 

6 T 

by Valenzuela. Ghia developed the backscatter using the Kirchhoff tangent- 
plane approximation, but real autocorrelations based on wave spectra were used 

Q 

rather than the usual mathematically convenient function. Jackson also used this 

method, but took into account the curvature, thus permitting extension of the 

9 10 

range of validity of the Kirchoff method. Semyonov, Barrick and Peake, 

11 

and Krishen used the Kirchoff method for part of the spectrum and added in 

the contribution, mostly important over a different angular range, of the small 

t 2 

perturbation theory; later Semyonov added a correction factor for Interrela- 

13 l4 

tions between the two parts of his solution. i?\mg and Chan s " developed a 
method for modifying the Kirchoff method in terms of the small perturbation 
solution for the undulating surface, thereby coupling the results from the two 
methods applicable to different parts of the wave spectrum and angular range. 

Regardless of the theory used, for incidence angles greater than about 
25° the primary factor governing ocean radar backscatter is the amplitude of the 
capillary portion of the sea spectrum. In the small perturbation theory, these ‘ 
components are readily related to the radar wavelength by the Bragg scattering 
condition that the radar return from adjacent crests of the most important 
ocean component adds in phase. This leads to the condition: 

A sin 0. - n X/2 

l 

where A = ocean wavelength (component of spectrum) 

X = radar wavelength 

0^ = angle of incidence (from local vertical). 
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The first order perturbation theory shows results only for n = 1, although 
intuitively one might expect higher order "Bragg resonances" also to be impor- 
tant. At an incident angle of 30° this means that the "Bragg wavelength" for 
the ocean is the same as the probing radar wavelength. Consequently, the Bragg 
wavelengths at that angle are, for the frequencies considered here , in the 
capillary region of the ocean wave spectrum. 

Measurements of capillary waves on the open sea are extremely difficult, 
for Che capillary waves usually are superimposed on other sea components many 
orders of magnitude larger. Consequently , our knowledge of capillary waves is 
based on theory and on measurements in wind-wave tanks . Recent work at HYU 15 * 16 
has combined the data from several investigators so that the variation in the 
capillary spectrum with wind speed can be determined. The data from three dif- 
ferent sources correlate well. 

The wave spectrum in the capillary region can be expressed in terms of 

boo = ik -° 5 *, 10 ~ 3) p 

where k is the wave number in cm ^ and 

D is a variable depending on wind speed . 

Experimental data indicate that S(k) continues to grow with increasing 
wind speed above about 3 m/sec. Below that speed, the capillary waves die 
away extremely rapidly. Thus, it is felt that the radar measurements will be 
invalid because the signals received will be dominated by returns from small 
patches where gusts are above the 3 m/sec. critical wind speed. 

Although preliminary analysis indicates an increase of backs catter with 
wind and wave development, none of the theories of electromagnetic backscatter 
presently account for breaking waves, sea foam, whitecaps and flying drops of 
spray . Under these conditions , sea surface is no longer a single valued func- 
tion of position and time. The development of an adequate theory to account 
for these effects is difficult. Good observations of sea return are rare for 
high winds although the available data show an increase in backscatter up to 
20 meters /sec. or so. 
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Measurements of Radar Backscatter 


Radar backscatter has been measured both in the laboratory and over the 
open ocean for a range of wind speed. Wright of the Waval Research Laboratory 
(WRL) has performed backscatter measurements at X-Band in a water tank with 
wind-driven waves as a function of wind speed and fetch. Preliminary results 
(as reported by Pierson and Stacy" ) indicate that the scattering coefficient 
c° at 60° incidence and vertical polarization increases with wind speed. 
Furthermore, Pierson claims that these results are quantitatively consistent 
with the growth of the capillary portion of the wave spectrum as a function of 
increasing wind speed. The range of wind speeds in Wright’s experiment cor- 
responds to ocean winds at 20 meter height of approximately 3m/s to 30 m/s. 

Perhaps the most extensive effort for measuring the radar sea return from 

17-20 

aircrafv was carried out by URL. Measurements were obtained as a function 

of polarization (vertical and horizontal}, incidence angle, and azimuth angle 

using four pulse radars operating at ^28, 1228, kk^5 , and 8910 MHz. Other 

21 22 

aircraft measurements * taken by WASA-JSC at 13.3 GHz used a fan-beam 

doppier radar technique to obtain a^. for incidence angles between ±60°. 

Wear nadir data were excluded. Recently measurements have also been obtained 

23 2k 

by WASA-LaRC, using a long pulse (beam limited) 13-9 GHz scatterometer * 

(AAFE RADSCAT) installed in a WASA-JSC aircraft. 

25 

An analysis of HRL and WASA-JSC data has been performed by Apel to 
determine the dependence of scattering coefficient on ocean surface wind speed. 
In an attempt to remove instrument biases in absolute cross section, he noimal- 
ized the value of cr° at a given incidence angle to the value of a° at 
6 i - 10°. A semilogarithmic presentation of the normalized differential 
scattering coefficient o°(35° )/o° (10° ) versus wind speed is given in figure 1. 
This type of data presentation has led to considerable controversy concerning 
whether or not the scattering coefficient saturates at the higher wind speeds. 
When these data are plotted on log-log paper (cross section in dB versus 
10 log^(wind speed}) as in figure 2, the result suggests a power law relation- 
ship such that 
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10 log <T° = 10 l,j A+ V(10 log U) 


0°(dE) = A(dB) + \)(10 log U) 


■where 


A - is a constant 
U - wind speed 

\> - is the wind speed power coefficient 

However, because of instrument biases and differences in operating frequency, 

20 

it is better to consider each instrument data set separately. Daley has 
curve-fitted power law response to NHL scattering coefficient versus wind 
speed data (wind speeds approximately 5 to 25 meters /second) for constant 
incidence angles. These power law curves for vertical and horizontal polari- 
zation are shown for upwind and downwind conditions in figures 3 and H respec- 
tively. The V’s for upwind and downwind observations for several incidence 
angles are given in the following table. 


— 

Freq. 

— 

Direction 



Polarization 

— 

Expo- 

nent 



Incidence Angle 

o 

O! 

■ ■ - 

25° 

30° 

U0° 

50° 

8.9 GHz 

Upwind 

W 

V 

0.20 

0.25 

0.37 


0.73 


Downwind 

W 

V 

0.20 

0.29 

0.36 

o.8o 

o.8o 


Upwind 

HH 

V 

0.00 

0.33 

0.58 

0.87 

1.03 


Downwind 

HH 

V 

0.00 

0.29 

0.52 

l.oU 

1.30 


Typical results from the NASA- JSC 13.3 GHz scatterometer (as analyzed by 
Bradley^ and Class sen^ of The University of Kansas) are shown in figure 5 
for the upwind case. The normalized differential scattering coefficient 
(a°(35 0 )/a°'{l0 o ) ) indicates a power law response of" wind speed to the l.hh 
power. The average exponents inferred from these data are given in the fol- 
lowing table. 
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Freq. 

Direction 

Polarization 

Exponent 

Incidence Angle, 0. 

25° 

35° 

13.3 GHz 

Upwind 

W 

\) 

1.12 

1.49 


Downwind 

W 

V 

1.15 

i.6o 


Crosswind 

W 

V 

1.00 

l.4o 


23 24 

Recently , preliminary results 5 have "been reduced from three aircraft 
flights of the AAFE RADSCAT. Upwind, downwind and crosswind observations of 
o and are presented versus wind speed for 0°, 10°, 20°, 25°, 30°, 40°, 

50°, and 55° in figures 6 and 7. The inferred wind speed exponents are given 
in the following table. 


Freq, 

Expo- 

nent 

Polarization 

Direction 

Incidence Angle, 0^ 

20° 

25° 

30° 

4o° 

50° 

13.9 C-Hz 

V 

W 

Upwind 

Downwind 

Crosswind 

1.0 
-99 
• 99 

1.53 

1.51 

l.5k 

1.9 

1.9 

1.9 

1.9 

1.89 

1.9 

1.9 

1.9 

1.9 

13.9 GHz 

V 

HH 

Upwind 

Downwind 

Crosswind 

.9k 
.9k 
• 76 

1.1*8 

1.1*8 

1.24 

1.9 

1.85 

1.69 

2.0 

1.98 

1.95 

2.0 

1.98 

1.95 


Previous observations indicated an azimuthal variation of cross section 
at the large incidence angles; i.e.: a wind direction dependence. To examine 

this effect the RADSCAT antenna was pointed to the nadir and then the aircraft 
flown in high banked circles. This maneuver simulated a conical scan of the 
ocean's surface. The results for three flights are shown in figure 8. Max- 
ima were obtained in the upwind and downwind directions with the upwind value 
being the greater. Minima were obtained for crosswind observations. The- 

shape of the curve is in qualitative agreement with Bragg scattering from the 































anisotropic capillary ■waves. These results imply that "both wind speed and 
direction may he inferred from multilook radar observations. Satellite mea- 
surements of scattering coefficient have been recently obtained during the 
NAS A- JSC Sky lab program. The S-193 RADSCAT experiment was a mechanically 

scanned pencil-beam 13*9 GHz scatterometer (s imi lar to AAFE RADSCAT). Fre- 
* „ 27 

limrnary results of ocean scattering coefficient versus wind speed are 
shown in figure 9 for horizontal polarization at 50° incidence angle. These 
data seem to verify the power law response to wind speed as do the measurements 
at other angles and polarizations. 

During several Skylab S-193 passes, simultaneous measurements were ob- 
tained with the AAFE aircraft instrument. The Skylab and AAFE RADSCAT had 
nearly identical antenna characteristics , but because of the large altitude 
difference the footprint for the satellite sensor was larger by greater than 
2 decades. Results from a pass on June 5, 19T3, in the Gulf of Mexico are 
shown in figure 10a, b, for vertical and horizontal polarizations respectively. 
The scattering coefficients for S-193 are approximately 2 . 5 dB greater than 
those for AAFE RADSCAT for incidence angles between nadir and 30°. Beyond 
30° the Skylab measurements fall off more rapidly than do the aircraft data 
which results in a crossing of the two curves around 1*0° incidence angle. 

While neither the AAFE or dkylab data are in final form the agreement is 
considered good. 


PERFORMANCE REQUIREMENTS 

During Phase A studies and in the early part of Phase B, the SeaSat study 
teams considered many candidate sensors to provide the User co mm u n ity with 
scientific and operational data in the areas of oceanography, meteorology, and 
geodesy. As a result of these studies, the microwave scatterometer was selec- 
ted as the prime wind vector sensor. Several configurations of each instru- 
ment were considered so that trade-offs could be made in designing the overall 
payload. The scatterometer options included fan beam antenna systems, scanning 
pencil beams and mixtures of the two. Decisions have now been made which nar- 
row the scatterometer options to a fixed fan beam antenna design. 
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Scatterometer specs have been developed for a prototype quasi-operational 
instrument to provide -wind speed and direction measurements (Mode I). A 
research mode is included to resolve uncertainties about the effects of sea 
state on wind sensing ability and to verify instrument design boundaries, 
such as optimum incidence angles and wind speed measurement range. In the 
development of these specs only wind speed measurement error due to instrument 
errors in scattering coefficient were considered. 

Frequency selection for the scatterometer has two conflicting directives. 

To minimize weather effects low frequencies (X-Band) are best while higher 

frequencies (Ku-Band) are preferred for maximum sensitivity to local wind 

speed. Skylab data taken at 13*9 GHz (Ku-Band) show favorable wind speed 

sensitivity which compares well with aircraft data also taken at 13.9 GHz. 

Less wind speed sensitivity is generally shown by HRL data taken at 8,9 GHz. 

Both of these points were discussed in detail in an earlier section. Recent 
28 

wave tank data confirm that Ku-Band should be more sensitive to surface 
winds than are lower band frequencies. Another factor concerning frequency 
selection is hardware development under AAFE RADSCAT, Skylab, and other pro- 
grams which also leads one to favor operation in the Ku-Band. 

Frequency allocation and RFI have been considered jointly with other 
onboard radars and application has been made for a scatterometer frequency 
between 14.5 and lU.9 GHz. The SeaSat altimeter frequency allocation request 
is for the 13. ^ to 13.9 GHz band. 


Mode I 

The scatterometer specs for Mode I were established from User require- 
ments (ref. l) of wind vector measurement range and accuracy, swath width, 
resolution cell size, and grid spacing (fig. ll). Orbit altitude, which is 
dictated by several factors, affects the scatterometer performance primarily 
by controlling swath width. The impact of other User requirements on the 
instrument specs will be seen in the discussion developed below. 
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The manner of transforming wind speed requirements to 0° accuracy (thus 

SCAT specs) can he seen hy examining the data in figure 12.*^ These da- 

ta are from the AAFE RAD SCAT aircraft instrument and correlate well with S-193 

Skylab data as shown in an earlier section. The allowable User wind speed error 

bar of ±2 m/sec or ±10 percent (whichever is greater) is projected on the a° 

curves for different incidence angles. By translating this error bar to the 0° 

axis* as shown, the corresponding instrument errur in o° is thus determined for 

each incidence angle and wind speed. These data are tabulated in fig. 13. At 

small incidence angles (less than 25° } the dependence of o° on wind speed is 

much less and the 0„ r exceeds the state of the art in measurement accuracy of 

M 2 

scattering cross section (±0.5 dB) 0 M /o° “ K = .12* ). In addition, for these 

m p 

incidence angles, 0° is also a function of sea state. For these reasons, earth 
incidence angles below about 25° are thereby excluded defining the inner edge 
of the measurement swath at 201 km from satellite subtract for an altitude of 
8o8 km (fig. l4), The outer edge of the main swath is set at about 55° earth 
incidence which is approximately the outer angle of present aircraft and S-193 
data. 

An additional high wind measurement zone (fig. 14) between 55° and 65° is 
specified since returned signals should be large enough to measure for high 
winds and it is expected that 0° will remain wind sensitive. The total 
swath width for high wind measurements is j6l km on each side of the satellite 
subtrack giving a total swath width of 1522 km, with an edge-to-edge sensitivity 
of 1924 km. 

Although the 0° sensitivity to wind speed variation is about the same 
for upwind, downwind, or crosswind conditions, the magnitude of 0° for cross- 
wind is less; hence, the crosswind design used here is a worst case design. 

In addition, since the instrument design is controlled by the expected 0° 
levels for low wind speeds, its performance as a wind sensor will generally 
exceed requirements for moderate and high wind speeds over the entire swath. 


“Vertical polarization 0° characteristics are used throughout this 
spec development whenever User requirements are involved directly. Horizontal 
polarization is also an instrument requirement. 

*2 

See link analysis section for definition of K . 

P 
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The scatterometer r s ability to measure wind speed does not cut off sharply* 
at low wind values but simply becomes less accurate. This is in keeping with 
User requirements since interest in wind speed begins to -wane in the 3 to 5 
m/sec zone; the scatterometer is specified to measure wind speeds at least as 
accurate as ±2 m/s ec down to 4 m/sec with decreased accuracy below that. The 
required values to accomplish this are shown in figure 15 s along with their 

respective o° values. 


The lower boundary of scatterometer dynamic range (figure 15) is based on 
measured a5t crosswind values for 4 m/sec (fig.. l6a). a° values are shown 
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in figure l6b for comparison. It should be reemphasized that the instrument 
will measure winds below 4 m/sec for the inner swath section but with reduced 
accuracy. The upper end of scatterometer dynamic range is based on a maximum 
upwind 0^ value of +10 dB determined by extrapolating curves in figure 6 to 
50 m/sec. For incidence angles below 15°, CF 0 ^ occurs at low wind speeds. 

Sines scatterometer calibration is achieved for 0° vs wind speed after 
orbiting the satellite, absolute a° measurements are not required. A fixed 
bias error in the absolute value of cr° would not affect wind speed accuracy 
because that bias would be calibrated out. Serious errors in a 0 would how- 
ever affect instrument measurement range and a ±2 dB accuracy has been selected 
to avoid this problem. The incremental system resolution of 0° needs to be 
±.25 dB so that variations in o° can be measured down to 10 percent without 
interference due to system limitations. This 10 percent standard deviation in 
0° is considered to be the lowest reasonable error for relative a° based on 
present satellite scatterometer designs capability. 


Since the ocean radar scattering coefficient is a function of wind 
direction as well as wind speed*, both a forward and an aft looking beam are 
required to obtain data at two azimuth angles for each resolution cell. An 
optimum implementation is forward and aft beam, each squinted 45° off the sub- 
satellite track to provide measurements which are separated in azimuth by 90°. 
In this configuration, the mean of the two measurements for a given resolution 
cell is nearly Independent of wind direction. The quadrant of wind direction 
can he determined from conventional meteorological data and satellite obtained 

*The dependence of a° on wind direction has been clearly demonstrated 
(figure 10 ) in the AAFE RADSCAT program and must be accounted for. 

Capillary* wave generation does cease abruptly at still lover wind speeds. 
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cloud mosaics and thereby improve the accuracy of the wind speed determina- 
tion. AAFE RADSCAT data have been processed in this manner to demonstrate 
feasibility (figure 8 and IT); however, there still are substantial amounts 
of experimental measurements required to define o' 0 wind direction dependence 
for a wide range of wind speeds and sea states. At present it is difficult 
to theoretically model this azimuth angle effect with any degree of confidence. 
A first effort has been made by Pierson (ref, 26) based on available aircraft 
data (figure l8a and b). Having determined wind speed from cr° values taken 
at orthogonal azimuth angles , individual forward and aft beam data can be used 
to determine wind direction. This is done by matching these two data points 
for a given resolution cell, fig. 19, to a previously determined experimental 
curve for that wind speed and incidence angle. This wind direction informa- 
tion can be used if necessary to iterate to an improved wind speed accuracy. 

A review of wind direction measurement errors using this two-azimuth 
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angle technique is given by Pierson . Generally, it shows the scatterometer 
can easily make ±20° accuracy measurements in wind direction for the higher 
wind speeds, but is less accurate for the lower wind speeds. There is no need 
for alarm, however, since certain mode options to be discussed in later sec- 
tions of this report give improved accuracy. 

In principle then, each resolution cell must be viewed from two azimuth 
angles to remove the wind direction effect and the best angular separation 
appears to be 90° in azimuth. The time between illumination of a given reso- 
lution cell by the forward and aft beams depends on the cells position along 
the fan beam illumination. Thus, the scatterometer must be designed properly 
in order to make the forward and aft beam cells cross. Each 50 km resolution 
cell must always have two footprints in it, thus, giving data at azimuth 

angles 90° apart as shown in fig. 19, and discussed in the subsequent design 
section. The doppler cell spacing along the fan beam is picked so that the 
cross ..track grid spacing for scat measurements are 100 km center-to-center 
and the along track grid spacing should be set at 100 km by the scat controller 
to comply with User requirements . In a high resolution mode the number of 
doppler cells along the fan beam could be increased to give 25 km resolution 

+ 

Private communication with W. J. Pierson, CUNY. 
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with essentially contiguous coverage. Tills is an option not presently planned 
for SeaSat A since it is more costly and User requirements do not specify the 
higher resolution. 

It should be pointed out, as shown in figure 11, that both the resolution 
cell size and grid spacing requirement impact the allowable scatterometer 
integration time and signal to noise requirements . 

Both and 0 ^ measurements are specifically required although pre- 

sent data indicate either one polarization or the other could satisfy User 
wind field measurement needs. In order to establish which polarization is 
best as a function of incidence angle and wind speed, both polarizations are 
included. Another aspect of dual polarization is that differential measure- 
ments may prove to be valuable for sensing ocean parameters. 

Polarisation purity of the wave incident on the Earth’s surface is required 
to avoid interference due to returns at the opposite polarization. This requires 
attention not only to spacecraft attitude and antenna alignment , but also to 
polarization purity as well. According to aircraft data, 0^. values can be 
as much as lOdB higher than 0 ° values for large incidence angles. In order 
to maintain the interferring vertically polarized signal 30 dB "below the hori- 
zontally polarized signal during the a° measurement, a polarization purity 

mi 

of 99 percent at the Earth’s surface was chosen along with a cross polarization 
antenna requirement of 20 dB. This combination assure", the 30 dB separation 
required between the two signals which limits that error source to about .25 dB 
(A a°/a° x 100 — 6 percent). 


Mode II 

The scatterometer specifications for the research mode have been estab- 
lished based primarily on instrument sponsor requirements. These requirements 
fall into two categories: One is concerned with studying parameters other 

than local winds (for example, sea state) which might affect the scatterometer 
performance and the other concerns data collection so that more definitive 
scatterometer design boundaries can be established (such as optimizing inci- 
dence angle and wind speed measurement range). Secondary requirements for 
Mode II includes the measurement of wind speed (not direction) along the 
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satellite sub track to supplement Mode X -wind speed measurements. Wind direc- 
tion would have to he derived from other sources such as synoptic wind fields 
developed from Mode I data taken concurrently or from ground truth data. 

Performance specifications of the scatterometer for Mode II are given in 
Fig. 20. The 0° measurement range for Mode II is extended "beyond Mode I 
specs to allow data taking from 0° to 65°. The fan beam is required to he ' 
oriented along the satellite subtrack so that (7° data can be taken (0° to 
65°) for a given patch of ocean. Highly instrumented test sites would be used 
whenever possible to give quality ground truth data. The actual angles are 
selected so that a smooth curve can be drawn to show the 0° dependence on 
angle for a given set of surface conditions. The steepness of the curve 
(ref. 29) warrants smaller angular steps near nadir. Scatterometer dynamic 
range and" cr° accuracy required for mode II was determined in the same manner 
as was used for Mode I. The 0° values listed near nadir are compatible 

0ISL3C 

with data taken with the Sky lab S-193 altimeter (ref. 30 ). Since Mode II 
required o° data for the full set of incidence angles from the same patch of 
ocean-, the only azimuth angles that are allowed are 0 and l80°; i.e., either a 
forward or aft antenna looking along the satellite subtrack. 

Polarization purity (99-7 percent) of the wave incident on the Earth's 
surface is specified for Mode II due to the cross polarisation measurement 
requirement. The high polarisation purity requirement was included to improve 
the accuracy for the cr° and a° measurements , where the return is lower 
than at higher incidence angles. There is no reliable data to indicate 

how much depolarization can occur at the ocean surface so the antenna cross 
polarization spec is die Gated by present antenna capability rather than by a 
rigid error analysis. It Is hoped that Skylab S-193 data will improve our 
knowledge in this area prior to the SeaSat launch but it is unlikely that such 
knowledge would warrant a major antenna development to improve the SeaSa+. A 
scatterometer polarisation purity. 

The geometric separation between measurements along the satellite sub- 
track is controlled by specific incidence angle selected and the timing of 
transmission from the Mode II antennas. If Mode II is being used to provide 
wind speed at 100 km grid points (for a specific incidence angle) thon timing 
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of the Mode'll transmissions must he adjusted accordingly. Data taken to pro- 
vide mnd speed should be restricted to incidence angles between approximately j 

25° and 55°. j 

| 

The resolution cell size for Mode II (50 km) assures compatibility with ] 

Users requirements when taking wind speed data and is small enough to allow . | 

ground truth data taken at the instrumented test sites to be representative j 

of the general conditions over the 50 km footprint (i.e,, no significant wind . 1 

I 

vector gradients across the footprint). : jj 

Other Modes 

Other operational modes are indicated in Figure 21. The actual modes of 
the scatterometer instrument will depend on the instrument design and it is 
not intended that all modes listed in figure 21 be considered as absolute 

requirements , nor to limit the potential modes; in fact, some of these modes \ 

may not be achievable. What is intended is that the instrument have some mode 

i 

flexibility beyond basic Modes I and II so that information derived from the 

early part of the SeaSat mission can be used to improve operation of the in- | 

st rumen t after launch of SeaSat A. 

The operational and scientific value of scatterometer data can be enhanced 
by judiciously using the instrument in various combinations of Modes I and II. 

For example, if only one polarization is used for Mode I, then portions of 
Mode II data could be taken in the time slot normally used for the other 
polarization. This still provides the User with wind speed and direction over 
the full design swath on both sides of the subtrack and rives wind speed along 
the satellite subtrack as well. 

Other combinations of Modes I and II have been selected (fig. 21 ) which 
will give flexibility to use of the instrument after launch. Either Modes. Ill 
or IV can be used after it is clear from early mission data which polariza- 
tion is best for operational use. This then leaves more instrument time to be 
given to data collection in the research mode (antenna 5* figure 1*0. Modes 
V through VIII allow improved a° measurement accuracy, with increased reso- 
lution cell size for one 500 km swath on either side of the satellite subtrack 
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(but not both). Mode IX and X allow maximum operational data t akin g for a.n 
five antennas and one polarization. 

A FAIT BEAM SCATTEROMETER DESIGN 
Mode I Analysis and Design 

To meet the specifications and design requirements discussed in the pre~ 
vious section a fan beam, design was chosen. This design incorporates four 
"stick" antennas for Mode I which produce a star-like illumination pattern on 
the surface of the earth. See Fig. l4. Backseatter data are sampled by each 
antenna with both vertical and horizontal polarization. With each polarization 
of each antenna following in turn, a scatterometer measurement is made from 
one antenna at a time. The data are obtained by transmitting a short pulse 
(5 msec) of CW signal then receiving, detecting, and integrating the return 
signal reflected from the ocean surface. Doppler filtering and range gating are 
used to divide the return signal along the fan beam into separate data cells at 
different incidence angles. An illumination pattern is produced which results 
in measurements over a wide range of incidence angles, with the forward looking 
and aft looking beams occurring at 90° azimuth with respect to each other. A 
half power heamwidth antenna of .5° x 25° was chosen as a compromise between 
antenna size and gain, swath coverage, and resolution capability. Figure 
illustrates the coverage geometry produced by the fan beam design for incidence 
angles of approximately 25° to 65° in Mode I. 

Figure 22 shows the side view geometry of one fan beam along with appli- 
cable equations defining the various geometric parameters, and figure 23 shows 
a plan view. Included in this figure are a few representative diagrams of 
resolution data, cells along one antenna beam illumination. Figure 2k shows 
an enlarged view of the geometry of a resolution cell. 

The instantaneous field of view (IFOV) or doppler cell is determined by 
the antenna heamwidth (narrow dimension), and doppler frequency line spacing 
(doppler bandwidth). The radar scattering coefficient (a 0 ) is determined 
from the radar return signal for each doppler cell averaged for a specified 
measurement period, (t^), resolution data cell on the surface of the 

21 


earth is then determined by the area generated by the doppler cell as a result 
of satellite motion during the measurement period. The equivalent resolution 
cell is defined by the square whose dimensions are the average of the dimen- 
sions of a rectangle within which the area generated by the moving doppler cell 
(IFOV) fits diagonally, as illustrated in figure 24. 

The orientation of the doppler lines within the beam illumination varies 
along the total beam as illustrated in figure 23. For the inner cells (low 
incidence angles) the doppler line is oriented approximately 45° to the central 
beam axis. For the outer cells (high incidence angles), the orientation is 
about 13° with respect to the beam central axis. In addition the width of the 
beam illumination increases as a function of slant range from the satellite 
(see figure 23). Figure 25 is a plot of a few constant doppler lines created 
by the satellite motion. Two beam illumination patterns are included in this 
figure to illustrate the orientation of the doppler lines with the antenna 
beam. In figure 24 the constant doppler lines forming a doppler cell are 
shown as straight parallel line approximations. In actuality, as seen in fig- 
ure 25 the constant doppler lines are approximately hyperbolas on the surface 
of the earth and are not exactly parallel-. Also it is assumed that the edges 
of the beam illumination are parallel inside a doppler cell, when actually 
they diverge slightly as pointed out previously. These approximations are 
made to simplify the calculation of the cell area (A c ), and the parameter L. 
These approximations should have negligible error effects on the overall 
system design. 

To establish a particular fan beam scatterometer instrument design. Certain 
parameters and constraints must be fixed. The parameters that are fixed include 
satellite altitude, resolution (r), spacing between successive scans of the 
same antenna beam (S ), number of antennas and polarizations to be sampled per 
scan, appro xim ate range of incidence angles (0.) to be sampled, transmitter 
frequency (f ), antenna half power beamwidths (£,$}, minimum and maximum 0° to 

"b 

be measured, and the maximum error allowable on the radar scattering coefficient 
(cJ°) measurement. In addition the center of each resolution cell of the aft 
antenna must overlap the data from the center of the corresponding resolution 
cell of the forward antenna. Using these parameters the remaining parameters 
which establish the resolution data cell geometry and the instrument character- 
istics can be determined. 
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Once the antenna L-s&mwidth is estahl: shed the choice of* doppler offset 
frequency of the IF and the corresponding handwidths must he determined. These 
frequencies along with the data sample period (t^) establish the resolution 
and geometric location of the doppler cells on the surface of the earth. To 
determine these frequencies the parameters 1 q (distance on earth surface 
from center to center of doppler cells along antenna "beam) and L (length of 
base of' doppler cell) must he calculated (see figures 23 and 2k ) . The follow- 


ing equations, derived in appendix 1, define these parameters, 

a 1 / 2 

^ vT 


i TT „ -1 

1 c 180* R COS 


/180 s 'C 


2 S s 


v t R $ sin <5 

, L = v5~r - R d, - ■ V' 


( 1 ) 

( 2 ) 


where : 

R = radius of the earth , km 

S - spacing along direction of satellite track between successive scans 
s 

of the same antenna, km 


r = equivalent resolution size, km 

v - satellite ground velocity, km/sec. 
g 

t = data measurement period, seconds 

p 1 

= slant range to center of cell, km 
$ - narrow 3 dB heamwidth of fan beam antenna, radians 

5 = angle between doppler line and line normal to satellite track, 

(doppler line rotation angle), degrees 

The dimension 1 is only a function of the antenna scan spacing (S ). This 
c s 

spacing of the resolution data cells must he picked to give cross track spac- 
ing equal to along track cell spacing. By choosing the proper earth angle (y^) 
to the first doppler cell the corresponding resolution data cells of the for- 
ward and aft antenna can he made to overlap at the center of the cell.* As 
noted in Appendix I, the earth angle y for any cell is a function of the 
instantaneous spacing (S c ) between corresponding forward and aft antenna 

doppler cells, and S is a function of the beam scan spacing (S_). Using 

c s 

*With earth rotation effects neglected 
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the folio-wing equation, -which is derived in Appendix 1, the proper value of 
the earth angle to the first doppler cell (y^) can he calculated. 



The variable 1 is a more complex parameter which varies in length along the 
beam illumination. Since the width of the beam illumination increases with 
incidence angle and the doppler line orientation changes, the value of L is 
made smaller for higher incidence angles in order to maintain the equivalent 
resolution cell constant along the beam illumination. If the parameter L is 
made a constant, then the equivalent resolution cell for higher incidence 
angles increases considerably. As seen from equation (2), L is a function of 
many parameters which include the equivalent resolution cell size, antenna 
beamwidth, slant range to the doppler cell, satellite ground velocity, data 
measurement period and doppler line rotation angle (see figure 24). Once 1 
and L for each cell are determined, the remaining geometric parameters, 
doppler frequencies and doppler cell bandwidths for each resolution data cell 
can be determined using the equations derived in Appendix 1. 


Measurement Timing 

There are two categories of timing sequences that must be considered: 
(l) the measurement-to-measurement timing and (2) the transmit /receive/ 
processing timing during each measurement. 


Measurement-to-Measurement Timing .- The time between the start of successive 

measurements from each beam scan is equal to S /v . Since there are four 

s g 

antenna beams in Mode I, each requiring two polarization measurements, then 
S /8v seconds is available per measurement period (t ). For the baseline 

s g P . , 

design, which is discussed later in more detail, S = 100 km, and v - 6.61 

s s 

km/sec. for an altitude of 8o8 km giving a t of 1.89 seconds. Using a 

t of 1.9 seconds, figure 26 shows one complete measurement sequence. Periodi- 
P 

cally this sequence would be interrupted for calibration measurements. Note 
in figure 26 that the beam sequence is 1, 4, 2, 3, so that data scans alternate 
between antennas on opposite sides of the satellite sub-track. This sequence 
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is the "best timing to allow center overlap of corresponding forward and aft 
data scans, and to allow the incidence angle of the first doppler cell to fall 
at approximately 25 °. 


Intra-Measurement Timing .- The following are the parameters of intra- 
measurement timing which must he determined for each doppler cell in the fan 
heam pattern. 

1. Transmit pulse length, 

2. Range gate "OR" time 

3. Range gate "OFF” time 

4. Noise gate "ON" time 

5. Noise gate "OFF" time 

6. Time between successive transmit pulses, T 

Figure 2J shows conceptually the timing for the transmit pulses, range gates 
for each cell, and noise gate times. The range gates are the periods when the 
signal plus noise from a particular doppler cell is filtered, square law 
detected, and integrated. The range gate time for each doppler cell is adjus- 
ted to match the two-way propagation time for each doppler cell. The noise 
gate is the period for filtering, detecting and integrating the mean noise 
level component with no signal present. The mean noise level is subtracted 
from the signal plus noise level measurement. This could be done by using a 
bi-polar integrator whose polarity is reversed during the noise integration 
period or by storing the noise level measurement for subtraction at a later 
time. To maximize the integration time, the transmit pulse length, 
should be picked as long as possible. However, the sum of T 

t 

transmit fall time, plus the i.f, relaxation time must be less than the two- 
way propagation time of the nearest point of the closest doppler cell. 


T , 
P 

plus the 


Once t has been chosen the range gate on and off times can be deter- 
P 

mined. The range gate on time is determined by the two-way range to the 
nearest point of a particular doppler cell. The range gate off time is 
determined by the sum of the two-way propagation time to the farthest point 
of that same doppler cell and the transmit pulse length The noise gate 

period is chosen to minimise the required signal-to-noise ratio (SNR) needed 
to obtain a specified measurement accuracy. As shown in Appendix 2, the 
minimum SNR occurs when the noise integration period is approximately twice 
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the sighal plus noise integration period. The noise on-and-off gate times are 
chosen to occur after all signal from a particular doppler cell has "been 
received. 

The time., T, "between successive transmit pulses is determined "by the 

transmit pulse length, the two-way range to the farthest point of the outer 

doppler cell, the length of the noise integration /period, and an additional 

small period of time to allow for data readout and pulse rise and fall times. 

The transmit puls s are repeated every T seconds throughout the measurement 

period, t . The total effective integration period, x_, T , as indicated in 
P 

figure 27 , is equal to the duty factor l-rf-j times the measurement period, t . 


Link Analysis 

To perform link calculations the following equations, which are derived 
in Appendix 2, are used. 


SNR. = 


R 


in KL B 
S c 


P (G/G ) 2 a°Le 2 X 2 L 

p — ° 

E “ r^ r * e 2 


(h) 

(5) 


SNR^ = signal noise ratio in the i.f. filter prior to the square law 
detector. 

Fp —true value of the radar return signal power, watts 
K = Boltzmann's constant 

T s - system noise temperature, degrees Kelvin 

B c = i.f . bandwidth of the filter preceding the square law detector. 


nominally set equal to the doppler spread over a cell, Hz 
- peak transmitter power, watts 
G = gain of antenna to center of doppler cell 
G q “ peak antenna gain 

cj> — narrow 3 d3 beamwidth of fan beam antenna, radians 
3 = wide 3 dB beamwidth of fan beam antenna, radians 
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0° - normalized radar scattering coefficient 

R c = slant range to center of doppler cell, Inn 

L - length of doppler cell along fan 1363111 illumination, tan 

£ = antenna efficiency 

X - wavelength of transmitted frequency, tan 

L g = miscellaneous system losses 

In equation ( 5 ) all parameters are essentially known except cj°. There- 
fore, the equation 

G° = C P R (6) 

can he written, where C = the value of all constants in equation (5). The 
true value of all the nonrandom parameters which make up C will not be 
known exactly, as well as the exact gains in the scatterometer receiver. Bias 
errors on the measurement of a° will occur as a result of not knowing these 
parameters exactly. Fluctuation or random errors will occur because the radar 
return signal from the surface of the earth is a Gaussian random process which 
is corrupted by system noise which is also a Gaussian random process. Both 
the bias and random errors will limit the absolute accuracy of the cr° 
measurement. However, the minimum change in a 0 measurement that can be 
resolved will be limited only by the random errors. The minimum resolution 
of a measurement is generally accepted as the standard deviation 0^ of the 
measurement , i . e . , 

a M - (Var a°) 1/2 (?) 

where Var 0° is the variance of the measurement of a°. Equation (6) can be 
written as , 

0° = C P R (8) 

A fk 

where a° is an estimate of 0° based on an estimate of the true value 

XV 

of return signal P_. If it is assumed that all parameters in equation (8) 
are essentially nonrandom except for P R , then is given by 


2T 


( 9 ) 


°M “ C ^ V 


1/2 


■where Var P„ is the variance of the measurement of P_. It is shown in 

n 

Appendix 2 that, 


K - 
P 


(Var P R ) 1/2 


( 10 ) 


•R 


where K is a normalised parameter which gives the standard deviation of the 
3? 

estimate of P R as a proportion of the true value of P fi . Combining equations 
(9), (10) and (6) yields. 


a.. ~ C K P~ - K a c 
M tj R p 


( 11 ) 


that is , the standard deviation 0,, of the measurement of 0° can he evalu- 

M 

r> 

ated directly hy determining the standard deviation of the measurement of P. 
The following equation, derived in Appendix 2, yields the parameter K 


R- 


■H?)' - (rt) (* * 2 “ * S H*kD 


( 12 ) 


where : 


k T B 

g G 

If/S = — ■ - — ■ = noise to signal ratio at the i.f. filter prior to square 


SR 

T R 

K. 


R 

law detection. 

= integration time of radar return signal plus noise, seconds 
- integration time for noise only, seconds 

= T R^ T SR 


Using equation ( 4 ) , ( 5 ) and (12), the expected standard deviation of the mea- 
surement of 0° can he determined for each doppler cell as a function of 0°. 
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Trade Off Analysis 


A trade off analysis was conducted to determine tile optimum choice of 
antenna beamwidth for a minimum transmitter power and highest equivalent 
resolution. This was done by chosing an antenna beamwidth and equivalent reso- 
lution and determining the transmitter power required to meet User requirements 
for swath width, minimum wind speed, and wind field measurement accuracy. A 
range of incidence angles of 25° to 55° was used in the analysis. Reasons for 
this choice of angles are given in the scatterometer performance requirements 
section. To cover this range of incidence angles a beam width of 25° in the 
wide dimension was required. Beamwidths in the narrow dimension of .35°, -5°, 
and 1° were used in the trade off analysis . Calculations were made for a 
range of equivalent resolutions of 10 to 60 km, with antenna scan spaeings of 
50 and 100 km, using four antennas and two polarizations. 

Table I lists some of the parameters and constants used in the analysis. 
The level of two-way system losses, L , was assumed equal to - 6.^ dB. This 
loss is made up of 3 dB losses in both transmitting and receiving and .H dB 
two-way loss through the atmosphere. The high transmitting and receiving 
losses are due to the required uses of rf switches and circulators, as seen 
in figure 38. The effective system noise temperature, T , referred to the 

£3 

input to the rf tunnel diode pre-amp, is equal tG. 1100 Kelvin. This value 

is determined using an assumed antenna brightness temperature of 200 K, a pre- 

amp noise figure of 6 dB, and a receiver loss of 3 dB. The value of T . is 

s 

calculated using the following equation: 


T 

s 



f 2 - (L r F‘~ 1} 
l r r 


(13) 


where , 

T = antenna brightness temperature, K 
& 

L = receiver loss, a factor larger than one 
K 

T q = reference temperature - 290 K 
F = tunnel diode pre-amplifier noise figure 

Combining and rearranging equations ^ and 5, the following equation 
results : 
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This equation -was used in the trade off analysis for calculating required 
transmitter power. The value of U/S used in this equation is determined 
from equation 12. Substituting in the value of .5 for K and solving for 
N/S, the following equation is obtained. 



Since = 2T gff> Xhis equation simplifies to, 



(15) 


(16) 


The terms B , R , and L, for a data cell are determined for each resolution 
c c 

considered in the analysis, using the equations of Appendix 1. The term 

g 

(G/G^) is determined from the antenna power pattern shown in figure 2 of 
Appendix 3, for the 25° beamwidth dimension. The values of the remaining 
parameters in equation (l4) are listed in table I. 

To minimize the required transmitter power the peak of the antenna beam 
must be directed toward the doppler cells with the higher incidence angles. 
This is done for two primary reasons; one, to compensate for the decrease in 
radar scattering coefficient, a°, at the higher incidence angles; and two, 
since the value of B c is smaller for the higher incidence angles, to pro- 
vide the increased SNR required to meet a .5 measurement accuracy. For the 
incidence angles considered in the trade off study (25° to 55°) s the antenna 
was pointed at an incidence angle of about yX ° , (0 - 43-5°) 5 so that the 
same transmitter power is required for both the 25° and 55° incidence angle 
measurements. This location was stermined by a few trial and error 
calculations. At this position the antenna gain is approximately - 9 dB down 
from the peak antenna for the 25° cell, and -.15 dB down for the 55° cell. 
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Using equation (l^) the required transmitter power -was calculated for 
various equivalent resolutions, antenna scan spacings, and antenna narrow 
dimension beamwidths. The results of these calculations are shown in figure 28. 

As seen from figure 28 the transmitter power required increases rather 
rapidly with decreasing resolution. This is due to the rapid decrease in the 
parameter L required to obtain the specified resolution, therefore reducing 
the bandwidth, B c , and increasing the SNR required to obtain a .5 measurement 
accuracy. For a fixed antenna beamwidth the parameter L is the only param- 
eter which can be reduced in order to obtain a smaller resolution. Since a 
narrow beamwidth dimension of the antenna produces a fixed illumination dimen- 
sion on the surface, there is a lower limit to resolution size. The limiting 
resolution for the 1°, ,5° s and .35° antenna beamwidths are approximately 35, 
15, and 10 km., respectively. 

As seen from figure 28, antenna beamwidths of .5° or smaller will be 
required to obtain resolutions of 50 km or smaller using reasonable trans- 
mitter powers. From antenna size considerations it appears that antenna beam- 
widths much smaller than .5° would be impractical (see Appendix 3). To meet 
the Users requirements for 50 km resolution an antenna beamwidth of 25° x .5° 
using a transmitter power of approximately 85 watts can be used to make U m/s 

wind speed measurements with a measurement accuracy (K ) of .5 over incidence 

P 

angles of 25 to 55 degrees at an altitude of 808 km. No system margin was 
assumed in this analysis. 

Baseline Scatterometer Fan Beam Design - Mode I 


After completing the trade-off study discussed in the previous section, 
a baseline design analysis was performed using the following design 
specification: 

Equivalent resolution low - 50 km; high - 25 km 

Antenna scan spacing low resolution - 100 km 

(along and cross track cell high resolution - 50 km 

spacing ) 


Swath width 


1000 Inn primary zone 
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Satellite altitude 
Wind speed range 
Accuracy- 

Range of a° 

1 standard deviation measure- 
ment error. K 
P 

Range of incidence angles 
Antenna beamwidth 
Number of antennas 
Number of polarisations 

(Non-rotating 

Low Resolution Design 


- 808 km 

- if m/s to if 8 m/s 

- ±2 m/s or 10 percent whichever is 
greater for incidence angles of 25 ° 
to 55° 

0 cLB to -28 cLB 

- 50 percent 

- - 25° to 65° 

- .5° x 25° 

- U 

- 2 (horizontal and vertical) 
earth assumed) 


With these specifications, calculations for all cell parameters were made. 
The value of 1 which provides proper overlap of the forward and aft data 
cells is 70.709 km. The two-way- propagation delay to the nearest cell is 
slightly under 6 msec., so a transmit on time as long as 5*7 msec, could, be 
used. However, in Mode XI operation, which will be discussed in a following 
section, the two-way propagation delay time to the 0° incidence angle doppler 
cell is approximately 5*3 msec. Therefore, for design simplicity, the transmit 
on time (t ) of 5 msec, was chosen for both Mode I and II. The integration 
time for noise only was chosen to he equal to twice the signal plus noise 
integration time. The following is a list of the various measurement times: 

t = 1.891 seconds 
P 

T -5 msec. 

P 

T = 26.8 msec. 

T-- t = 352.8 msec 
= 705*6 msec 

t = 1.0 msec, 
m 

Duty Factor = 19 percent 


To provide center overlap of the corresponding aft and forward data cells 
the location of the first doppler cell is required to be placed at an inci- 
dence angle of 2k. 1 degrees (see equation 3). 
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Using the equations in Appendix 1 all geometric parameters for each data 
cell were calculated for the low resolution "baseline design. These parameters 
are listed in table II. A diagram of the data cells on the surface of the 
earth from one forward and aft antenna beam is shown in figure 29. The range 
gate on and off times for each cell is summarised in table III. 


As seen from figure 29, a total of 15 contiguous data cells are produced 
within a fan beam illumination. To provide the 100 km along- and cross-track 
grid spacing between each data cell, every other data cell is used. This is 
illustrated in figure 30, where a more detailed diagram of the data cell 
illumination on the surface of the earth is presented. Measurements from 
additional surface area can be obtained by including the remaining T doppler 
cells (shown by dotted outlines in figure 30). Using these additional cells, 
a 50 km by 50 km spacing of data cells is provided. System performance is not 
affected by including these measurements. However, a slight cost increase 
will be incurred for additional doppler filters, detectors, and integrators. 


To determine the measurement accuracy capability of the baseline design 

for various wind speeds , a link calculation for each data cell was performed. 

To do this a calculation of the noise to signal ratio (H/S) for each cell 

was made using equations (k) and (5). Substituting these values of N/S into 

equation (12) the measurement accuracy, K , was determined for various wind 

P* 

speeds for the incidence angle corresponding to each data cell. 


The parameter a° for each wind speed (crosswind, vertical polarization) 
and incidence angle are obtained from figure 6c or figure 15 * 

The parameters B c and L are taken from table II, for the corresponding 
data cell, Hote that the value of L for cells 11 through 15 was made equal 
to the value of L at cell number 10, This was done to maintain a reasonable 
value of doppler bandwidth, B^, since the value of L becomes very small for 
these outer cells, when trying to maintain a 50 km resolution. The resolution 
of these outer cells will be bigger than 50 km. For example, cell number 11 
will have a resolution of approximately 53 km, and cell number 15 a resolution 
of TO km. 

The term (G/G^) is usually determined from the antenna radiation pattern 
such as that shown in figure 2 of Appendix 3. However, for these calculations 
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the approximate equations, 

o /0 - 0 , 

(G/e o ) 2 = -2 4 z - g" -i 


3 


dB 


(17) 


where , 


0p = the nadir angle from the satellite to the direction of the peak, 
antenna gain, and 

0^ = the nadir angle from the satellite to doppler cell, n, was used to 

determine this parameter for each data cell. The peak gain of the antenna is 

pointed slightly beyond cell number 7 (6 ~ 42°). This is done, as explained 

P 

in the trade-off study, to provide maximum gain for the outer doppler cells, 
where the scattering coefficient cr° is the lowest, and thereby minimise the 
required transmitter power. Due to limited spacecraft dc power and TWTA 
availability, a maximum transmitter power of 125 watts has been chosen for the 
baseline scatterometer design. The values of the remaining constants -In equa- 
tions (4) and (5) are listed as follows: 


K = 1.38 x 10" 23 , joules/ K 

T = 1100° 
s 

<f> = 8.75 x 10 radians, (.5°) 

3 = .43 6 radians, (25°) 

£ - • 5 

A - 2.157 x 10 3 km, ( f^ - 13.9 GHz) must be updated for SeaSat 
frequency allocation (approximately l4.6 GHz) 

L - -6.4 dB, (-3 cLB each for receiver and transmitter losses and -.4 dB 

5 

for two-way atmospheric losses) 


The results of the link calculations for the lowest a° design values 
corresponding to a 4 m/s wind speed (crosswind, vertical polarization) are 
summarized in table IV. The actual instrument accuracies, resulting from the 
link analysis, for data cells whose incidence angles correspond approximately 
to those required (figure l6) are summarized in table V, for various wind 
speeds. A review of tables IV and V will show that a measurement accuracy of 
.5 or better was only obtained for cells 3 through 9 (6^ = 83° to 53°) for a 
wind speed of 4 m/sec., thus cutting down the swath width at which 4 m/sec. 
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wind speed measurements can "be made. However, at slightly higher wind speeds 
(6 m/s) the specified measurement accuracies can he extended to 0^ = 24° to 
56° (500 km swath). At a wind speed of 3 m/s an accuracy measurement of 50 
percent can only he obtained for some of the middle data cells (0^ = 33° to 
48°). The location of the peak antenna gain can he changed to greatly improve 
the performance of the low incidence angle data cells at a slight sacrifice 
in the performance of the higher incidence angles. By doing this the 4 m/sec. 
measurement swath width can he extended. If the peak antenna gain is placed 
at a 0^ of 4l°, the measurement accuracies on data cells 1 and 2 are improved 
sufficiently to obtain 4 m/sec. measurement capability from cells 1 through 9 
(450 km swath width). A minimum windspeed measurement capability of 5 m/sec., 
and 12 m/sec. is obtained at cells 10 and 15 respectively. 

A comparison of tables IV and V with figure l6 shows that some margin 
exists in the measurement accuracy over those required, for the middle inci- 
dence angles at low wind speeds, and for the low and high incidence angles at 
higher wind speeds. It is generally desirable to provide a certain amount of 
margin to cover changes in the value of various parameters, such as transmitter 
power, transmitter and receiver losses, and receiver gains. Figure 31 is a 
plot which shows the amount of signal margin obtained as a function of wind 
speed for the incidence angles shown in table IV. 

Figure 32 is a plot showing how the measurement accuracy varies as a 
function of the scattering coefficient, a 0 , for the J incidence angles listed 
in table V. Similar curves which lie in the region between cell number 1 and 
number ? would result for the remaining eight data cells. This figure shows 
that the normalised measurement error, K^, is a monotonically decreasing func- 
tion of a°. This is to say, it also decreases as signal-to-noise ratio (SNR) 
increases, since a° is proportional to SNR as seen from equations (4) and 
(5). This relationship is shown in figure 33 where a 0 is plotted as a func- 
tion of SNR for three data cells. Also shown plotted is the normalized 

measurement error, K , as a function of SNB. Note that for SNR < 0 dB, K 

P 

decreases proportional to the increase in SNR. For SNR » 0 dB, no significant 
decrease in can he obtained by increasing SNR. This is due to the inher- 

ent Gaussian fluctuation of the radar return signals and the limited number 
of measurement samples taken of the ocean surface. The equivalent number 
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of sample measurements taken can be shown to be the product of the doppler 

bandwidth and integration or sample time, that is, the product B T a „. low 

this limit affects the ultimate measurement accp^^cy can be seen from equation 

(12), where K approaches a value of 
P 

The ultimate accuracy varies, since eac 
width. For cell number 1 the ultimate accuracy is approximately 1 percent, 
and for cell number 15 this value is slightly over 3 percent. 

The measurement accuracy at any SNR can be improved at the expense of 
resolution by averaging the measurement of one or more adjacent data cells. 

This is the same as increasing the number of sample measurements or the 
B c T SN Pr° duct > i'k can *> e seen effectively increases the resolution size. 

The measurement accuracy will improve by the square root of the number of data 
cells averaged. If data cells along the antenna beam are averaged the resolu- 
tion size will effectively increase by a factor equal to the number of cells 
averaged. However, if successive scans of the same antenna are made so that 
adjacent cells are contiguous, then averaging along track and/or cross, track 
data cells will increase the resolution size by only a small percentage deter- 
mined mainly by the sample period, t^. This is more easily seen by a review 
of figures 23 and 24. 

A review of equation (4) and ( 5 ) will show, after a few sample calcula- 

3 

tions , that the term B R /L is essentially a constant for all data cells 

^ c 3 

(see table IV). This is due to the increase in the value of R being 

o 

directly offset by the corresponding decrease in the value of B c and L. 
Therefore, the only variable between data cells which affects the SHR received 
for a particular a° is the antenna gain function (G/G^) , since the remain- 
ing parameters in equation (5) are constants. Note also that for a particular 

SHR the value of K varies between data cells as a function only of the 

P 

square root of the ratios of their doppler bandwidths , B c « 

High Resolution Design 

A high resolution design analysis was performed in an attempt to obtain 
25 km resolution capability, at higher wind speed conditions. This was done 
by dividing the 50 km resolution doppler cells into two 25 km doppler cells. 


B T 
c l SN 


' 


, as SNR approaches infinity. 


h cell has a different donuler band- 
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so that 50 km resolution measurements could still be obtained by averaging 
measurements from the two 25 km floppier cells. If the value of the parameter 
It is reduced in half from Its value at 50 km resolution, received SHE will 
be the same since the bandwidth, B c? will also be reduced in half. However, 
the ® c T gjj product will also be reduced in half, resulting in the measurement 
error of one 25 km cell being increased by the VfF". This increase, however, 
is offset by the \^ 2 ~ improvement resulting from averaging two 25 km data 
cells as discussed In the previous section and Appendix II. Thus, no degrada- 
tion results In meeting the original 50 km resolution measurement accuracy 
requirements. 

A review of the equation defining the parameter L will show that It 
must he reduced by more than a factor of two in order to reduce the equivalent 
resolution from 50 to 25 km. In fact, for the inner cells, L must be reduced 
by a factor slightly greater than 3- For the middle and outer cells, L 
becomes very small because of the effect of the doppler line orientation. 
Because of this as noted previously in the 50 km baseline design, the equiva- 
lent resolutions of cells 11 through 15 were made greater than 50 km. For 
these outer cells, at best, L can be reduced by a factor of k; however, the 
resolution of these cells will still be about 35 to 53 km. Therefore, little 
will be gained by dividing the outer cells into two higher resolution cells. 
For the middle cells, L can be reduced by a factor of 3 or k, resulting in 
equivalent resolutions of 26 to 33 km. The measurement error of one 25 km 
cell relative to a 50 km cell will increase by /T, where w is the factor 
by which L is reduced from i’js 50 km value. Since an improvement of only 

is obtained by averaging two 25 km cell measurements (to obtain a 50 km 
measurement), additional transmitter power or slightly higher 0° (higher wind 
speed) , equal to /w/2 , will be required to obtain the same measurement 

accuracy as that obtained by using a single 50 km resolution data cell. 

With this approach a high resolution (25 km) baseline design was con- 
ducted by dividing each of the first ten 50 km cells into two 25 km cells. 

For this design a closer grid spacing of 50 km was used. For this spacing it 
is advisable to either scan the antennas at a higher rate, which would reduce 
the integration time, or to limit the number of antennas or polarizations 
sc ann ed. For the higher scan rate the decreased integration time will degrade 
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the measurement accuracy, limiting the data to higher wind speeds. In an 
actual system design both approaches could be incorporated, since it would 
involve only a control function change in the timing to provide either mode 
of operation. For the high resolution baseline design a four-antenna one- 
polarization approach was used, which provides a 1000 km swath width. Figure 
34 illustrates a sample of the scan pattern for this baseline design. As seen 
from this figure, a total of 20 contiguous data cells are produced within a 
fan beam illumination for incidence angles from 25° to 55°. As in the ease of 
the 50 km baseline design, every other 25 km data cell is used to provide the 
50 km along- and cross-track grid spacing (figure 3*0. 

The results of the calculations for 25 km resolution are shown in Tables 
VI and VII, for each incidence angle. Minimum a° with a measurement accuracy 
of .5, and the approximate corresponding windspeeds are shown. Also shown in 
these tables are the results obtained by combining two 25 km cell measurements 
to obtain a 50 km cell measurement. As seen from the tables, wind speed mea- 
surements of 4 m/sec can he made at design accuracies for incidence angles 
ranging from about 34 to 4l degrees for 25 km resolution and about 33 to 51 
degrees for 50 km resolution. At the lower and higher incidence angles the 
full swath width requirement is met only if windspeed exceeds 3 m/s for the 
25 km resolution case, and 6 m/s for 50 km resolution case. 

Mode II Scatterometer Design 

Mode II is the experimental operating mode in which a fan beam antenna 
is directed along the satellite subtrack to measure 0° incidence angles 
of 0 to 65 degrees for the same ocean area. The general specifications and 
design requirements for Mode II operation are discussed in a previous section. 
The basic area of illumination for the Mode II antenna is illustrated in 
figure 14. Designs for this mode are described in the following sections. 

Doppler Cell Design 

Definition of the dependence of radar scattering coefficient on 0^ 
require a° measurements every 5° of incidence angle from 0 to 65 degrees 
plus a measurement at 2.5°. In addition measurement.' at these incidence 
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angles must cover the same specified surface area. There are two possible 
design approaches which can be used to meet these requirements. The first is 
to choose doppler cells at the exact incidence angles required. In this 
approach the center of the cells would not be equidistant on the surface of 
the earth’, therefore, measurements at the required incidence angles would 
require variable timing to cover the same surface area. The second approach 
is to choose doppler cells whose center-to-center surface spacing is some 
fixed increments or a multiple of this increment. The measurement time inter- 
val would then be equal to this spacing or its multiple divided by the satel- 
lite ground velocity. In this approach each incidence angle measurement 
would cover the same surface area; however, the incidence angles would not 
occur exactly at the required values. 

Figure 35 illustrates the geometry of the surface illumination for the 
two options and one polarization. Also included in the figure is an enlarged 
view of the resolution data cell. The resolution data cell is determined 
partly by the doppler bandwidth, which establishes the parameter L, and the 
smear distance generated during the measurement period. Any combination of 
doppler cell length {L) and measurement period (t^) can be chosen to form a 
specified resolution length. It is shown in Appendix 2 that the optimum 
choice of these parameters, which maximizes the produetjis made when 

the smear length and the doppler cell, length are made equal. A 50 km resolu- 
tion length has been specified; therefore, the doppler cell length and the 
smear length have been set equal to 25 km. With a satellite ground speed of 
6.6l km/s , the measurement period, t , of 3.782 seconds is obtained. 

In figure 35 option I shows the fixed increments of incidence angle 
approach with variable surface spacing of the data cells. Also note the over- 
lap of the 0° and 5° cells by the 2.5° cell. Option II shows the fixed sur- 
face spacing approach, using increments of 50 km along with their correspond- 
ing incidence angle. The data cells are contiguous for the lower incidence 
angles. At the higher incidence angles only those cells are included whose 
incidence angles are closest to the required values. 

A measurement made at the other polarization would he shifted oh the 
surface of the earth by 25 km. Therefore, the data cells of the two 
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polarization measurements would only overlap the same surface area by 25 km 
(50 percent). The two additional eross-polarization measurements would not 
overlap the first polarization measurement and only part of the second polari- 
zation measurement. A total of 125 km would be covered by the four polariza- 
tion measurements. To provide better overlapping coverage of the four 
polarization measurements . the measurement period, t p . could fce reduoed and 
the doppler cell length (L) increased. However, as shown in Appendix 2, the 

B T ow product would he greatly reduced resulting in a less accurate measurement, 
c orJ 


Measurement Timing 

Intr a-Me asur ement Timing . - In order to simplify the electronics, the 
timing philosophy is kept the same as Mode I, The transmit pulse length in 
Mode II is therefore kept at 5 ms. Also the pulse repetition period was made 
equal to 26.8 msec., the same as Mode I. For option I, the range gate on-times 
are shown in table VIII. The range gate off-times are 5 msec later. Option II 
timings are approximately the same. 

Me asurement-to-Me asur ement Timing . - To produce overlap of the same sur- 
face area by each incidence angle measurement, the time interval between 
measurements must be property adjusted. 

Figure 36 illustrates the measurement sequence required for both options 
at H and V polarization for each incidence angle measurement to cover the same 
surface area. As seen from the figure for both option I and II it requires 
over 200 seconds and a surface distance of over 1300 km for all incidence 
angles to measure the same surface area. At each incidence angle a 3*78 
second measurement is made for each of two polarizations. Sufficient time is 
not available to make cross-polarization measurements that cover the same 
surface area. Cell number 15, (Q i = 65°) will be the first cell to cover a 
specified surface area, followed in sequence by the remaining lower incidence 
angles. The times shown between each measurement are the times required to 
allow the spacecraft to traverse the appropriate distance so that the measure- 
ment at the next incidence angle will occur at the specified surface area. 

The times are shown for a 808 km altitude spacecraft with a ground speed of 
6,6l km/s. Note that in option I the 2.5° incidence angle (cell number 2) 
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is not included, since for two polarisation measurements, sufficient time is 
not available between cells 1 and 3 for the 2.5° measurement to cover the 
specified area. With one polarization measurement however, cell number 2 
could cover the same specified area. It must be noted however, that during 
each 3-7$ second measurement period, measurements for all incidence angles are 
being obtained including an incidence angle of 2.5°, even though they may not 
be at the same specified area. 

The major difference between option I and option II is the sampling 

interval used between measurements. In option II the measurement interval is 

the time it takes to move 50 km on the surface, or multiples thereof. As suchj 

completely new measurements could be made each 50 km {t = 7 . 56 s,), thus 

P 

allowing each incidence angle to contiguously cover the same surface of the 
earth. In option I the time intervals are not at a fixed rate,' therefore, 
all incidence angles will not necessarily fall on the same surface of the 
earth except for the surface area measured by cell number 15* 

During the 200 plus seconds that Mode II data are being obtained. Mode I 
will be inoperative, since both these modes will be using common receiver 
front ends and transmitter. Certain mode combinations are desirable, such 
as those shown on figure 21. For example, some Mode I data can be obtained 
in the surface area adjacent to that of Mode II. This can be accomplished if 
only one polarization is used in both modes. This is illustrated in figure 
37, where a combined measurement sequence for Mode I and Mode II using one 
polarization is shown. Mode II measurement timing is arranged to allow each 
incidence angle to measure the same surface area starting with cell 15 * 

During the off-time between each Mode II measurement period, four Mode I 
measurements are made, one for each antenna. The measurement timing for 
Mode I is adjusted to allow the same surface overlap of the front and back 
antennas for most of the measurement periods. Since the off -time interval 
of Mode II at the lower incidence angles is small, only two Mode I antennas 
can be sampled during those time periods (see figure 37). 

Link Calculation .- A link calculation was performed for option I. 

Equations h, 5, 15, and 16 were used to perform these calculations. Prior to 
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r 

i 

these calculations the geometric parameters for each doppler cell were first 
determined using the equations of Appendix 1. The spacecraft angle, 0, slant j 

J 

range, R , and doppler bandwidth, B , Obtained from these calculations are 
c c 

required in the link calculations. The parameter L is the same for all data 

cells and is equal to 25 km. The integration period for signal plus noise is 

705.6 ms., and for noise only is 1.411 sec. The remaining parameters in the * j 

i 

link equations are the same as those in Mode I except the antenna beamwidth, 

B. 

Since a relatively wide S/C angle, 0, of about 54° must be covered for i 

Mode II measurements, a wider beamwidth antenna than used in Mode I is j 

! 

required. The same antenna design as used in Mode I can be used with one 
dimension being changed to increase the beamwidth B, as discussed in 
Appendix, 3. A beamwidth of ,5° x 40° was chosen for the Mode II design. The 
40° beamwidth was chosen as a compromise on coverage, antenna gain loss and 
maximum beamwidth obtainable using a modified Mode I antenna design 

(Appendix 3). The beam center was pointed at S/C angle, 0 , of 42.1°. This i 

o ** i 

angle was picked so that the required minimum a measurement of 0 dB can 

be made at an incidence angle of 0°. 

Using equation 15, the N/S required to obtain a measurement accuracy j 

(K ) of .5 was determined for each data cell- Then using equations 4 and 5, 

P 

the O value which would provide this N/S value was calculated. Table VIII 

t 

lists the geometric parameters and the results of the link calculations. As 
seen from the table the minimum o° obtained at incidence angles above 50° 
did not meet the required minimums listed in Fig. 20. Since these calculations 
were made with a transmitter power equal to 125 watts, very little margin is 
available for some of the incidence angles which did meet the required 
minimums. The only alternative to obtaining the required minimum a° 
measurements, with adequate margin, is to use two .5° x 25° beamwidth ' 

antennas . 

Earth Rotational Effects | 

In the fan beam design presented in the previous sections the effects j 

j 

of. the earth’s rotation were not included. This was done to simplify the p 

i. 

initial design analysis. However, in the final system design the effects 
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of earth rotation on the doppler frequencies and bandwidths, location 
of data cells on the earth’s surface, and the measurement overlap requirements 
of the forward and aft antenna cells will have to be taken into 
consideration. A brief discussion of these effects is presented in this 
section. A more detailed analysis with applicable equations is presented in 
Appendix 1. 

At first it would appear that the tangential motion (.465 krn/s at the 
equator) of the earth's rotation would have a small effect on the fan beam 
design. This is not the case. Two effects, which are major consequences, 
occur on the system design: one is the change in the doppler velocity 

created by a component of the earth’s tangential velocity, and the other is 
the movement of the surface area measured by the forward antenna data cell 
during the time required for thg spacecraft to travel the surface distance 
(S^) separating the forward and aft antenna data cells. Each of these 
problems is discussed in the following paragraphs . 

The tangential component of the earth's rotational velocity must be 
combined with the satellite velocity component to establish the doppler filter 
frequencies and bandwidths. This component, although small in comparison to 
the satellite velocity and easily accounted for in doppler velocity 
calculations does not remain constant in amplitude or orientation with 
respect to the satellite velocity throughout the satellite orbit. For this 
reason the doppler filter frequencies and bandwidths for any particular 
incidence angle will vary throughout the orbit. Since the orbit is 
retrograde, the satellite velocity vector at the maximum latitude points 
(72° N and S) is 180° away from the earth's tangential velocity vector. 
Throughout the remaining portion of the orbit the satellite velocity vector 
(as well as antenna beams) rotates +72° from this position. The maximum 
rotation of 72° occurs at the equator. The earth's tangential velocity 
decreases with increasing latitude, from its maximum value at the equator, 
as a function of the cosine of the latitude angle. 

The total doppler velocity^ on all antenna beams is approximately equal 


^Doppler velocity in this discussion refers to the component in the plane 
of the antenna beam, parallel to the earth's surface. This velocity must 
he multiplied by sin0 to obtain the doppler velocity to a particular data 
cell along the beam (see figure 22) . 
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at the maximum latitudes of 72°. During the remaining portion of the orbit 
the doppler velocity is either increased or decreased from this value as a 
function of the changing orientation of the antenna beams with respect to 

the earth's velocity vector. During the ascending portion of the orbit the 

magnitude of the doppler velocity is increased on beams 1 and 3, and decreased 

on beams 2 and 4. During the descending portion of the orbit the situation 

is reversed. In both cases the magnitude of the doppler velocities on beams 1 

and 3 are always approximately equal. The same holds true for beams 2 and 4. 

The slight inequality that does occur is caused by the latitude difference 

between these beams at any point in the orbit. At the equator the latitude 

differences have a negligible effect on the doppler velocities. At the 

maximum spacecraft latitude (72°) the difference in latitude between beams 1 

and 3 and between beams 2 and 4 causes about a 1.7°/ difference in their 

o 

doppler velocities at a 65 incidence angle and about a .4 / difference at 
a 24° incidence angle. 

The doppler velocities on each beam change about 6°/ o as the spacecraft 
moves from the maximum latitude to the equator. Since this change is 
positive in one crossing direction and negative in the other crossing 
direction, a total change in doppler velocity of about 12°/ occurs on each 
beam throughout one orbit. Since only one bank of fixed doppler filters 
will be provided in the receiver processor, the effect of the variable 
doppler velocity will be to shift the location of the data cell corresponding 
to each doppler filter in and out along each beam. The shift will thus cause 
an expanding and compressing of the range of incidence angles throughout 
the orbit. For example, the incidence angle for cell number 1 will vary 
from approximately 22*7° to 25.8°, for cell number 10 from approximately 51° 
to 61 , and for cell number 15 from 60 to JP . 

One possible solution to alleviate the wide variation in incidence 
angles caused by the varying doppler velocity is to rotate the spacecraft 
approximately +3*3° about the yaw axis during a single orbit. It is shown 
in Appendix 1 that by slight continuous yawing of the antenna b^ams the 
doppler velocity can be held essentially constant throughout the orbit* The 
effect of this yawing is to compensate for the variable earth velocity vector 
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by varying the component of satellite velocity (function of cos 4* * » see figure 
22) which contributes to the total doppler velocity. Instead of a continuous 
yawing, a few step changes in yaw during specific portions of the orbit 
could result in significant improvement. The rotation of the spacecraft 
impacts the design of the satellite bus, and as such may be costly to 
incorporate. However, the scatter ometer can be utilised without rotation 
capability with a minimum impact on User requirements. 

The capability of providing forward and aft antenna data cell overlap 
at the same surface of the earth is greatly affected by the earth's rotation. 
At the equator, the surface of the earth at data cell number 1 of the forward 
antenna will have moved approximately 31 km by the time cell number 1 of the 
aft antenna reaches the corresponding location. Thus only a small corner 
of the equivalent resolution square of the two data cells will overlap. At 
the maximum incidence angle (65°) the surface measured by the front antenna 
data cell will have moved approximately 130 km before the aft antenna data 
cell reaches the corresponding location. The remaining data cell measurements 
of the front antenna will move by some amount between these extreme values. 

At the higher latitudes, of course, the amount of shift of the front antenna 
data measurements is correspondingly smaller due to the decrease in the 
earth's tangential velocity. (See Appendix 1 for a more detailed discussion 
of the amount and direction of data cell movement caused by the earth's 
rotation throughout the satellite orbit). 

Thus it is seen that the ideal situation of data cell overlap depicted 
in figure 30 is greatly altered by the rotational motion of the earth. This 
movement of the forward antenna data cells considerably compromises the 
ability to make orthogonal measurements to determine surface wind direction 
and to remove the effects of wind direction in making a wind speed measure- 
ment. A study of figure 30 will show that when the data cell measurements 
of the front antenna are shifted by the amount indicated in Appendix 1, 
they lie either in a blank area between data cells or they fall in an area 
measured by the aft antenna at a higher incidence angle - or a combination 
of both. A possible solution to providing orthogonal measurements of the 
same surface of the earth from the front and back antennas is to include 
all contiguous data cells along each antenna beam and make the along— track 
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spacing (S ) equal to 50 km. By doing this the blank spaces shown in 
s 

figure 30 will be filled in, therefore providing contiguous coverage of 
50 km squares within which scattering coefficient measurements are made. Thus 
two orthogonal measurement cells can be found in which common overlap area 
does exist, although the incidence angles at which these measurements are 
made are not exactly the same. (Maximum difference in incidence angle that 
will occur is approximately 4.5°). The sample time, t , will have to be 
reduced in half to provide 50 km along-track spacing. The measurement 
integration time will also be reduced in half which will result in a 
measurement accuracy degradation by /2. The minimum wind speed that can be 
measured will be slightly higher (about 1 m/s) than the baseline design 
capability. 

Mode II will be affected in a similar manner, in trying to obtain 
measurements at the same surface of the earth at all incidence angles. 

However, if the test sites are located at a high latitude (60° or higher) , 
the 65° incidence angle measurement will only move by about 40 to 50 km 
until the 0° incidence angle measurement is made. The higher incidence 
angle measurements will fall closer to the 65° measurement since there is less 
delay between those measurements. 

Data reduction and processing tasks will be significantly increased 
due to earth’s rotational effects. Precise orbital and earth-local parameters 
are required so that correct detennination of the location and incidence 
angles of the resoluticn cells can be made for each measurement, and so 
proper combining of forward and aft antenna cell measurements are assured. 

In summary then, the following changes to the baseline design can be 
made to accommodate the effects of the earth’s rotation: 

1. Provide contiguous data cells for each antenna beam and set the 
timing to provide along-track data cell spacing of 50 km. 

2. Locate sites for Mode II measurements at latitudes of 60° or greater, 
or at lower latitudes in groups of 5 incidence angles on successive orbits. 

BASELINE HARDWARE SYSTEM DESCRIPTION 
System Block Diagram 

The particular design considered here is based largely on systems used in 
other scatterometers , such as AAZE RADSCAT and Skylab S-193. Because of this, 
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some sections contain more detail than would otherwise he available. This 
design is intended only as an example. 

A block diagram of the SeaSat Scatterometer is shown in Figure 38. This 
system utilises an integrated frequency synthesizer to generate all of the 
necessary rf signals at UKF and microwave frequencies needed by the exciter, 
transmitter, mixers and scatterometer processor. The synthesizer consists of 
a highly stable crystal controlled oscillator from which all required frequen- 
cies are derived thus assuring coherency throughout the sensor. The low level 
microwave output of the synthesizer is amplitude gated (on and off) by a pin 
diode type switch to insure high isolation (typically _> 60dB) between the "ON" 
and "OFF" states. The pulsed rf is amplified in the final output stage which 
consists of a TWT and integrated power supply. The 125W peak power output 
pulse is ’directed to the stick antenna array through a. series of waveguide 
components consisting of electromechanical waveguide switches, circulators and 
1X5 waveguide diode switch. Each antenna has associated with it a switchable 
polarization circulator which is commanded in a preset fashion by the timing 
and control electronics. 

Received rf energy from the proper antenna is directed through the chosen 
path of the 5 -way switch and sent through the T/R circulator and receiver pro- 
tection circulator to the low noise preamplifier, band pass filter and first 
conversion mixer. In the mixer a microwave signal derived from the stable 
oscillator in the synthesizer is added to the return signal to produce a 
coherent IF signal which is then fed to the scatterometer processor. The Scat- 
terometer processor consists of IF amplifiers, amplitude=control range gates, 
and processing channels. The processing channels for the data consist of a 
bank of band-pass doppler filters, detectors, dc amplifiers and signal inte- 
grators. The timing and control for the processor are provided by the control 
unit. 

The function of the control unit is to accept spacecraft commands and 
power, to generate the precise timing and control logic needed by the scattero- 
meter to form rf pulses and operate the processor. In addition it accepts the 
output scatterometer data and instrument telemetry parameters and formats them 
so they are compatible with the spacecraft data system. 
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Scatterometer Transmitter 


Integrated Frequency Synthesizer .- Using an ultra stable crystal oscilla- 
tor the Integrated frequency synthesizer generates all required rf signals at UHF, 
VHF and microwave frequencies for the scatterometer. Stabilities of better than 1 
part in 10° over a few milliseconds and 1 part in 10 per day can he obtained 
with proportionally controlled heated crystal oscillators. The oscillator 
frequency chosen depends on the final output frequency selected as well as the 
desired 1st and 2nd rf frequencies, but would typically be 50 MHz or below 
and employ a 3rd overtone crystal. The basic oscillator output is then buf- 
fered, amplified, and used to drive a harmonic generator which produces inte- 
ger multiples of this frequency. Using stripline or combline filter techniques 
the required UHF/VHF signals for down conversion are then obtained. Another 
portion of the crystal oscillator’s output is used to drive a phased locked 
sclid state multiplier and produce the proper microwave signal for the output. 

Fin Diode Switch - This component utilizes a pair of pin diodes shunted 
across the waveguide structure which are driven into conduction or back biased 
into the "OFF” state, typical performance achievable in such a switch is 
> 35 dB OH/OFF ratio, < 2 dB insertion loss in OH state, <_ 10 nsec rise and 
fall time and <2 nsec rf pulse delay. These performance specifications insure 
high fidelity pulse with ultra sharp rise and fall times to guarantee proper 
scatterometer performance . 

Output TWTA .- The output TWTA amplifies the pulsed rf input from the pin 
diode switch and has an output of 125W minimum saturated power. Based on a 
TWT design derived from a Hughes tube (294-H) to be flown on the Japanese 

\ 

broadcast Satellite, the TWTA can operate over a 200 MHz bandwidth with greater 
than 40 percent efficiency including the power supply. Typical gain is h2 dB 
minimum at saturation^ and advanced cathode design should permit operation for 
over 3 years with little or no performance degradation. The TWTA requires 
spacecraft power input, commands to turn on the heaters and high voltage plus 
a pulse modulation input command. It provides telemetry outputs of critical 
temperatures, voltages and currents ■go the scatterometer housekeeping telemetry 
data stream. Although reliability figures for space qualified TWT’s are 
greater than 90 percent, overall amplifier system reliability can be further 
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enhanced "by parallel operation of two or more tubes which provides redundancy 
or can he used to increase total output power if operated simultaneously. 
However, in this simultaneous mode the relative phase of the TWTA output must 
he maintained within 25 ° (ref. 31 ) in order for the signals to add properly. 
Under present technology, there is no convincing evidence that TWTA r s can he 
operated together with this phase stability. 

Scatterometer Receiver 

Receiver "Front End” .- The scatterometer receiver front end consists of 
the 5-way waveguide diode switch which directs energy from, the chosen antenna 
to the T/R and receiver protection circulator. Both of these circulators 
together provide about 60 dB of isolation for the receiver when the pulse is 
being transmitted. In the receive state they should exhibit less than 0.5 dB 
insertion loss. The low noise device will be a Tunnel Diode Amplifier (TDA). It 
has a gain of 30 dB with a 6 dB noise figure. Minimizing losses from the TDA 
to the antenna should keep the total system noise temperature to about 1100 K 
(assuming total loss of 3 dB). 

The output of the TDA is combined with a low level microwave signal derived 
from the synthesizer in a low noise, 3ow conversion loss diode mixer to produce 
the 1st IF signal. 

Scatterometer Processor — The scatterometer processor must amplify, detect 
and integrate the samples of the backscattered energy and provide a measure of 
the total received power to the spacecraft data system. To do this, systems 
of this nature generally amplify the 1st IF signal, then mix once again with 
a signal derived from the synthesizer to bring the information down to a con- 
venient 2nd IF frequency. At this frequency the range gating and doppler 
filtering are accomplished which achieve the desired spatial resolution on the 
ground. After filtering, the selected element is detected, amplified and 
integrated for a specific period of time. It is this integrated value, which 
represents an est imat e of the received power, which is A/D converted and put in 
the data stream. 

The scatterometer processor is controlled through logic and switching 
commands provided by the digital control electronics to insure the proper 
sequence of range gating and integration. 
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Control Unit 

The function of this unit is to act as interface between the spacecraft 
and the scatterometer providing all necessary input voltages and accepting 
output data and housekeeping telemetry. A few of its functions are: generation 
of pulses of proper Width and PKF for scatterometer transmitter, generation 
of proper sequence of commands for switching antennas and antenna polarization, 
generation of proper commands for operation of scatterometer processor includ- 
ing range gates and integration commands , distribution of raw and processed 
voltages to scatterometer components, formating, A/D converting, and multi- 
plexing of input data and telemetry, and delivery of the scatterometer data 
stream to the spacecraft data handling system. 

Antenna System 

A total of five separate antennas are used in V a scatterometer design. 
Each antenna can operate in one of two linear polarizav r on modes. The four 
antennas used in Mode I have an antenna beamwidth of .5 ^ 25 degrees. The 
fifth antenna has a beamwidth of „5 x 40 degrees and is use? in Mode II only. 
Various antenna designs can be utilized to produce the required beamwidths. 

The design requires a smooth uniform antenna pattern in both planes , with low 
side lobe levels and high polarization purity. For a more detydl discussion, 
of an antenna design which produces these qualities see Appendix 3. 

CONCLUDING REMARKS 

A historical review of microwave scatterometer aircraft and satellite 
measurements has been given to show that the instrument can be used to sense 
ocean winds on a global basis. Instrument performance characteristics have 
been developed based on User requirements for a wind sensor on the SeaSat-A 
satellite. User requirements include a wind speed measurement range of 4-50 
m/sec ±2 m/sec for a swath width of 1000 km, a resolution cell size of 50 km, 
and a measurement grid spacing of 100 km. 

A specific instrument design was developed that satisfies the User require- 
ments and specifications with minor exceptions. This instrument uses 4 


50 


antenna arrays, each 'with a .5° x 25° half power heamwidth, producing a j6l km 
measurement swath on each side of the satellite. Wind speed measurements are 
made over a range of incidence angles from 24° to 66°. A mini mum wind speed 
measurement of 4 m/s can he made for incidence angles of 24° to 53° (450 km 
swath on each side) and a 5 m/s measurement for incidence angles of 24° to 56 ° 
(500 km swath on each side). Higher wind speeds (7 to 12 m/sec. ) are required 
for satisfactory measurement at the outer incidence angles. 

A research mode instrument design using a fifth antenna (heamwidth = ,5° x 
40°) was developed to make o° measurements along the satellite sub-track for 
incidence angles of 0° to 65°. However, the minimum a 0 measurement require- 
ments for this mode can not he completely satisfied for the higher incidence 
angles (50°to65°). The minimum a° that can he measured at an incidence 
angle of 45° is -24.5 dB, which is only 1/2 dB from the required levels. 

Changes to the baseline design can be made to accommodate the effects 
of Earth's rotation as follows: 

1. Provide contiguous data cells for each antenna beam and set the 
timing to provide along-track data cell spacing of 50 km. 

2. Locate test sites for Mode II measurements at latitudes of 60° 
or greater, or obtain measurements in groups of 5 incidence angles 

on successive orbits at a lower latitude. 
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Table I - Trade off analysis parameters and constants. 


a° = - 10.8 dB for 0^^ = 25°, and - 2?. 5 dB for 0 = 55° 

Altitude “ 808 Ion 

Measurement accuracy, K = .5 

Transmit pulse width, r = 5*5 msec. 

P 

Pulse repetition period, T = 26. ? msec. 

Integration period for signal plus noise, x^ - 195 msec, for 50 Inn spacing 
and 389 msec, for 100 km spacing 

Integration period for noise only, x^ = 2x giJ 

_5 

Transmitted frequency, f^ — 13*9 GHz (X = 2.15T x 10 km) 

Antenna efficiency, e = .5 
Antenna beamwidth, 0 - 25° 

Antenna beamwidth, <f> = . 35°j *5° and 1.0° 
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Table II - Baseline Geometric Parameters - MODE I 


£23 

1 

Y, cLeg 

6, deg 


1 

mu 

6, deg 

L, km 


■ 

318.26 

2.85 

21.30 

24.15 

875.76 

177377 

5. Cl 

52.75 

25002.8 

2 

388,97 

3.49 

25.36 

28.86 

907.34 

209199 

6.86 

51.78 

22055.1 

3 

459.68 

4.12 

29.H 

33-24 

943,84 

237611 

8-95 

50.55 


■ 

530.39 

4.76 

32.54 

37.31 

984.70 

262710 

11.18 

49.03 

16306.1 

5 

601.10 

5.39 

35.66 

4i.o6 

1029.39 

284713 

13.48 

47.22 

13725.7 

6 

671.81 

6.03 

38.48 

44.52 

1077.45 

303901 

15.76 

45.06 

11401.0 

■ 

142.52 

6.67 

41.03 

47.70 

1128.43 

320582 

17.99 

42.54 

9348.0 

8 

813.23 

7.30 

' 

43.32 

50.63 

1181.94 

335063 

20.13 

39.61 

7553.3 

9 

883.95 

7.94 

45.39 

53-33 

1237.66 

347630 

22.14 

36.24 


10 

954.66 

8.57 

47.24 

55.62 

1295*29 

358544 

24.03 

32.37 

4655.2 

11 

1025.37 

9.21 

48.91 

58.12 

1354.58 

368033 

mm 

32.37 

4051.1 

12 

1096.08 

9.84 

50.41 

60.25 

1415.32 

376299 

27.38 

32.37 

3532.0 

13 

1166.79 

10.48 

51.75 

62.23 

1477.31 

383512 

28.87 

32.37 

3085.6 

14 

1237.50 

11.11 

52.97 

64.08 

1540.41 

389820 

30*22 

32.37 

2701.3 

15 

1308.21 

11.75 

54.06 

65.81 

1604.48 

395348 

31.47 

32.37 

PMi 




Alt, - 808 km. 

r = 50 km. 

1 = 70-709 km 

S = 100 km. 
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f t = 13.9 GHz 



















































































































Table III Range Gate Timing - MODE I 


Cell Humber 

B 

2 

3 

4 

5 

6 

7 

8 

9 

R. G. ON (ms) 

5.-8 

6.1 

6.3 

6.6 

6. 9 

7.2 

B 

7-9 

8.3 

R. G. OFF (ms) 

10.8 

11.1 

11.3 


11.9 

12.2 

12.5 

12.9 

13.3 

E. G. WIDTH (ms) 

5 

5 

5 

5 

5 

5 

5 

5 

5 


10 

11 

12 

13 

8.6 

9.0 

B 

9.9 

13.6 

l4.o 

Ik.h 

14.9 

5 

5 

5 

5 


Ik 

10.3 

15.3 
5 


15 


10. 7 


15. T 
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Table XV Baseline Link Analysis and Accuracy Parameters For 4 ra/s Wind Speed - MODE X 


CELL 

Wo. 

BnlP 

L 

KT B 
s c 

B c T SN 

(G/G q )^ dB 

8 

O 

t> 

S/N, dfi 

It , 8/0 
P 

■ 

3.18X10 11 

3.8xi0“ l6 

8820.23 

-16.46 

-16 

-17.27 

70.4 

2 

3.18X10 11 

3.3X10" 16 

7780.37 

-10.63 

-21 

-16.44 

62.1 

3 

3.18X10 11 

2.9 X 10 -1 ^ 

67^2.36 

-6.38 

-22.5 

-13.69 

35.9 

n 

3.17X10 11 

2.5 x 10“ l6 

5752.31 

-3.44 

-24 

-12.24 

28.1 

5 

• 

3 . 17 x 10^ 1 

■j 

2.1x10 

4842.0 

-1.54 

-24.2 

-10.54 

20.2 

6 

3.16X10 11 

l.7*lcT l6 

1*021.93 

- .476 

-26 

-11.26 

27.1 

7 

3.16X10 11 

l.lxio" 16 

3297.89 

- .036 

-26.3 

-11.11 

29.0 

8 

3.15X10 11 

l.lxlO" 16 

2664.58 

- .058 

-26.8 

-11.62 

37*7 

9 

3.14X10 11 

9.1X10" 17 

2116.02 

- .441 

-27 

-12.19 

45.9 

10 

3.13X10 11 

7.ixio~ 1 ^ 

1642.21 

-1.05 

-28 

-13.79 

74.3 

11 

3.11X10 11 

6.lxiO _1T 

l429.ll 

-1.83 

-28 

-14.55 

94.5 

12 

3.09X10 11 

5 . 4xio“ iT 

1245.98 

-2.72 

-28 

-15-40 

123 

13 

3.07X10 11 


1088.52 

-3.65 

-28 

-16.31 

l6l 

l4 

3.05X10 11 

i+.ixio -17 

952.94 

-4.62 

-28 

-17.25 

213 

15 

3.02X10 11 

3.6xio“ 1T 

835.94 

-5.59 

-28 

-18.17 

28l 

. 



BW 


T 

s 

■wind 


= 125 watts 
= 42° 

= 1.891 seconds 

- 26.803 msec. 

- 5 msec. 

= 35^.77 msec. 

= 705-53 msec. 

= . 5° x 25° 

= 1100 It 
speed = 4m/ s 

















TABLE V FAN BEAM SCATTEROMETER ACCURACIES - MODE I 



0 = k2° 

P 

P,p = 125 watt s 

ant . Bw = . 5° x 25° 

t = 1.891 sec 
P 

T =26.8 msec 


T SN = 3 ^2.8 rasec 

T„ = 705.5 msec 
N 

r = 50 km 

S = 100 km 
s 

altitude = 808 km 


-i4»« 





















Table VI-Geometric Parameters and Minimum a° and Wind speeds for 25 km Resolution Design (Cells 1 throu gh 5) - MODE I 


50km. 
Cell Wo. 

1 

2 

3 

4 

5 

s s = 50km 

Alt. = 808 km 

P T = 125 v 

t = 1.891s 
P 

T SW = 352.8 ms 

25km 

Cell 

lower 

upper 



lower 

upper 

lower 

upper 

lower 

upper 

lower 

upper 

Rg, km 

300.5 

336.0 

371.2 

406.7 

441.9 

477.4 

512.6 

548,2 

583.3 

618.9 

6, deg 

20.2 


24.7 

26.3 

28.2 

30.0 

31.7 

33.4 

34.9 

36.4 

0i, deg 

22.9 

25.4 

27-7 

30.0 

32.2 

34.3 

36.3 

38.3 

40.1 

41.9 

S 3 km 

G 

868.7 

883.2 

899.0 

916.1 

934.2 

953.7 

974.0 

995-6 

1017*8 

1041.2 

L, km 

17-3 

17.3 

16.29 

16.29 

15.02 

15.02 

13.46 

13.46 

11.58 

11.58 

3 c> Hz 

8404.5 

7994.6 

7136.4 

6740.7 

5856.7 

5502.2 

4625-6 

4328.6 

3482.9 

3250.5 


-10.3 

-13.6 

-l6.2 

-18.6 

-20.4 

-22 

-23.1 

-24.2 

-24.6 

-25.1 

IC 

P 

m 

.5 

■ 5 

.5 

.5 

.5 

• 5 

.5 

.5 

• 5 

w/s,m/s 

9 

6 

. 

6 

5.5 

5 

4 

4 

3.5 

3.5 

3.5 

o°. , dB 
min 

-13 

-18.7 

-22*7 



-25.2 

-26.5 


.914 

.44 

.87 

.51 

.84 

• 58 

.79 

.63 

.76 

.67 


.5 

.'5 

.5 

.5 

.5 

v/s ,m/s 

6.5 

5 

. 

4 

3 
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TABLE VII - Geometric Parameters and Minimum <3° and Wind speeds for 25 km Resolution Design (Cells 5 through 10) 
MODE I 


50 km 

6 

7 

8 

9 

1 

0 

25 krn 
Cell 

lower 

upper 

lower 

upper 

lower 

upper 

lower 

upper 

lower 

upper 

R s km 
g 

65^.9 

688.7 

726.6 

758.5 

798.4 

828.1 

870.4 

897.5 

942.5 

966.8 

0, deg 

37*8 

39.1 

40.5 

41.6 

42.9 

43,8 

45 

45.8 

46.9 

47.5 

0.j deg 
i 

43*7 

45.3 

47 

48.4 

50 

51.2 

52.8 

53.8 

55.4 

56.2 

R , km 
o’ 

1065.7 

1089.4 

1116.7 

1140.3 

1165.5 

1193.5 

1226.8 

1248.6 

1285.3 

1305.4 

L, km 

11.27 

11.27 

10.64 

10.64 

9.9 

9-9 

9.1 


8.1 

8,1 


2947.6 

2756.7 

2413.3 

2263.9 

1944.5 

1832 

1547.2 

1465.7 

1193 

1137.1 

c £in dB 

-25.5 

-25.7 

-25.6 

-25.6 

-25.8 

-25-0 

-24,4 

-24,1 

-23.2 

-22.9 

K 

T3 

• 5 

.5 

• 5 

.5 

.5 

.5 

.5 

.5 

.5 

■ 5 

w/s, m/s 

3-5 

4 

4 

HI 

5 

5 

5.5 

6 

7 

7.5 

•in" 

-27 *2 

-27.2 

-27.0 

— 

-25.9 

-2! 

4.7 


.725 

.694 

.704 

.715 

— 

.647 

.776 

.695 

.745 

.691 

• 745 


.5 

• 5 

.5 

.5 

.5 

w/s, ra/s 

3 

3.5 

4 

5 

6 












































































Table VIII Geometric Parameters, Link Parameters and Minimum a° Measurements for Mode II (Option I) 


m 

B 

R , km 
S 

0, deg. 

G i ,deg. 

R -Inn 
c 3 

On Gate, 
ms 

f I)> hz 

B ,hz 
c 

KT B 
s^c 
wants 

— 

VsB 

(G/G o ) 2 

dB 

S/fl 

dB 

a° 

dB 

a 

0.0 

0.0 

0.0 

808 

5.39 

0.0 

21367 

}.2xi0“ 16 

15076 

-26.7 

-17 

0 

2 

31.31 

2.21 

2.5 

808.7 

5.40 

26740 

21312 

3.2X10 -16 

15033 

-23.9 

-16.9 

-2.7 

3 

62.71 

4.43 

5 

811 

5.40 

53428 

21150 

3.2X10" 16 

" 4923 

-21.3 

-16.9 

-5.2 

a 

126.2 

8.86 

10 

819 

5.45 

106446 

20511 

1 £ 

3.1x10 

14472 

-16.6 

-16.9 

-10 

5 

191.2 

13.28 

15 

833 

5.54 

156657 

19481 

3xio" 16 

13746 

-12.5 

-16.8 

_ 

-13.9 

6 

259.2 

17.67 

20 

854 

5.67 

209661 

18111 

2.8xl0~ ie 

12779 

- 9 

-l6.6 

-17.2 

a 

330.6 

22.03 

25 

881 

5.84 

259067 

16467 

-16 

2 . 5*10 

11619 

- 6 

-I6.3 

-19.9 

8 

l+06,8 

26.34 

30 

916 

6.07 

306498 

14627 

2.2x10"^ 

10321 

- 3.7 

-l6.1 

-22 

9 

489.5 

30,6 

35 

961 

6. 36 

351604 

12674 

1.9><Kf ie 

8943 

- 2 

-15.8 

-23.4 

10 

580.4 

34.78 

40 

1016 

6.72 

394031 

10693 

i.6xio" lf 

7545 

- .8 

-15 A 

-24.2 

11 

682 

38.87 

45 

1085 

7.18 

433458 

8764 

1.3xi0 -lt 

6184 

.16 

-15 

-24.5 

12 

797.4 

42.83 

50 

1170 

7.74 

469586 

6958 

1.1X10' 1£ 

4910 

- .01 

-14.5 

-24.1 

13 

930.7 

46,63 

55 

1276 

8.44 

502141 

5331 


3762 

- .31 

-13.8 

-23.2 

14 

1087.3 

50.23 

60 

l4o8 

9-32 

530875 

3924 

6xio -1 ^ 

n 

- 1.0 

-13.2 

-21.9 

15 

1274.3 

53.55 

65 

1574 

10.42 

555569 

2755 


1944 

- 1.9 



-12.4 

-20.2 


0 = 42.1° 

P 

t = 3.782 sec. 
P 

t SH = * 53 ms 

= l.Ull sec. 
Btf ** . 5°xl+0° 

s ■ - 5 

L = 25 km 

Remaining 
Parameters same 
as Mode I 









































UPWIND DATA 
■ MISSION 88 

• MISSION 119 13.3 GHz 
▲ MISSION 156 

♦ NRL (8.9 GHz) 


Wind speed, knots 

Normalized scattering coefficient data from JSC and NRL 
(log-linear plot). 





<r at a = i(r 



Wind speed, knots 

Figure 2.- Normalized scatlering coefficient data from JSC and NRL 
(log -log plot). 















'Wind speed, knots 


Figure 5.- Upwind JSC normalized scattering coefficient vs. wind speed. 



















Wind azimuth angle, degrees 

Figure 8. - AAFE Radscat Scatterometer wind direction sensitivity. 




Scattering coefficient 



Figure 9.- Sky lab S-193 Scat ter ometer response to 
wind at sea (from congressional presentation 
by Dr. R. K. Moore, Feb. 20 } 1974). 




Figure 10.- Comparison of AAFE Radscat with Sliylab S-i.93 data. 

June 5, 1973 pass, Gulf of Mexico, 13 knot winds, upwind. 
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Figure 11. DEVELOPMENT OF SCATTEROMETER SPECIFICATIONS 
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Figure 12.- Wind speed 
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{crosswind AAFE Fiadscat data). 






* These values are determined from extrapolation of curves beyond 


the 3 m/sec to 26 m/sec data range. 

Figure 13.- Required o° v accuracy for scatterometer wind speed accuracy 
of ± 2 m/sec or ± 10 percent. 
































High and Low Wind Swath 


TOP VIEW 



Figure 14.- Scatterometer fan beam geometry. 
















• Operating Frequency Bond 14,5 - 14.9 GHz 

• Scattering Cross Section 

• c° normalized Standard Deviation 



• Range of o° Values (See Note 5) 


v. incaa^/ice 
\An^Ie 

ParameteiN^^ 

23° 


40° 

50° 

55° 

6 o° 
Hote 1 

6 b° 
Note .! 

a° Max dH 

*7 

+6 

+5 

+2 

0 

0 

0 

o° Min dp 

-19 

-24 

-27 

-30 

-31 

-31 

-31 


• Maximum Mas Error + 2 dB 

• Resolution Requirement 0.2 dB 

• Antenna Requirements (See Note 5) 

• S;.de Lobe Level -20 dB 

• Azimuth Angles - Two required 90° apart. Data should be timed to causa the 

two resolution cell centers to overlap. 

• Earth incidence Angle 25° to 65° (See Note 1) 
fi Polarization Options - w and HII 

» Polarization Purity - 99 percent (At Earth Surface) 

® Cross Polarization level - 20 dB integrated value across the half power beam relative 
to integrated value of main beam primary polarization. 

• Geometric Separation Between i ‘ J Measurements 

• Along Track 100 km 

• Cross Track 100 km 

• Resolution Cell Size 50 km 

(including integration smear) 


Note 1, The minitmai measurable a 0 value at >0' and tb 0 is desirable but not 
mandatory. 

Note 2. These values are determined from exi rnpoiation of o°vv ■ ■u'ves beyond 
the present data range (3 to 26 m/nec. ) 


Note 3, The range of values is based on RAPSCAT aircraft data (upwind , downwind and 
crosswind) and has a 3 dB margin at both ends of the o° range. 


«**§* 


Note 4, Values are from figure 16a for a crosswind of 4 tn/sec. 

Mote 5, Signal return due to side lobes must be 30 dn below return signal from 
main beam for that some doppler band. 


FIGURE 15 BCATTEIOMETER PERFORMANCE SPECIFICATION FOR KOBE I 
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Incidence angle, degrees 

(b) I-IorizontaL polarization. 

Figure 16.- Concluded. 
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Figure 18.- 


Wind direction, degrees 
(a) Oyy vs. wind azimuth angle. 

Scattering coefficient deoendence on wind azimuth angle at 8 
(from reference 26). 
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(a) Earth coverage geometry and scatterometer doppler cells 



Figure 19.- Orthogonal doppler cell crossings for forward and aft antenna beams (Mode I). 


• Operating Frequency Same as Mode I 

• Scattering Grogs Section Measurement Ability 

• c° normalized standard deviation - optimize by proper choice of antenna 
beamwidth J3 and 6_ (Tube power fixed by Mode X requirements). 
Beamwidth angle <J> p must be same as Mode I to minimize antenna 
development cost. 

• Range of a° values (See Note 1) 


'~''^_Frequency 
^'"-'^^Angle 
Par amet er~-'^( 

— 

0 

10 

20 

25 

30 

40 

50 

55 

60 

65 

a °max » 

+20 

+9 

+5 

+ 7 

+6 

+3 

+2 

0 

0 

0 

cr° . , dB 

min 

+4 

+2 

-13 

-19 

~ 2 h 

-27 

-30 

1 

H 

-31 

-31 


• Maximum o° bias error h_ 2 dB 

• o° resolution requirement 0.2 dB 

• Antenna Requirements 

• Azimuth angle cP or l8o 1 ' relative to satellite subtrack vector 

• Earth incidence angle 0 to 6‘i 0 

• Side lobe level -20 dB relative to peak gain 

• Polarization options o vv , H , ^, lv , o v]I 

• Polarization purity 99 percent 

• Cross polarisation level -20 dB 

(Integrated value over the half power beam relative to integrated 
value of main bean, primary polarization) 

• Resolution Cell Size 30 km (including integration smear) 

• Geometric Separation Between o° Measurements 

• Mode II (Common resolution cell o° measurement for each incidence angle. ) 

• Mode III, IX, and X 100 km along track. 

Note 1. The range of values is based on RADSCAT aircraft data (upwind, downwind, 
and crosswind) and has a 3 dB margin at both ends of the a° range. 

Note 2. Signal return due to integrated effects of side lobes must be 20 dB 
below return signal from main beam for that same doppler band. 


FIGURE 20 SCATTEROMETER PERFORMANCE SPECIFICATION MODE II 
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Mode 

Identification 

Antenna 

Polarization 

Swath 

Width 


mEmmmSm 

Antenna no. 
5 

Comments 

I 

o o 

°W’ °HH 

— 

1000 Km 

Mode for talcing operational data 

II 


0 o O 0 

a VV> ^HH’ °VH’ a HV 

NA 

Research mode for evaluating instrument limits 

in 

°VV 

Any two pol . 

1000 Km 

Gives user req'd wind speed and direction plus 
wind speed along the satellite subtrack 

XV 

°HH 

Any two pol. 

1000 Km 

Gives user req r d wind speed and direction plus 
wind speed along the satellite subtrack 

V 

Oyy antennas 
1 and 2 only 

— 

50u Km 

Improves a 0 measurement accuracy for V , 
operations’ data poi 

VI 

Oyy antennas 
3 and 4 only 

- 

500 Km 

Improves a 0 measurement accuracy for V , 
operational data 

vn 

°HH an<:ennas 
1 and 2 only 

— 

500 Km 

Improves cr° measurement accuracy for H , 
operational data poi 

vm 

ff HH an * :<!Vnas 
3 and 4 only 


500 Km 

Improves a 0 measurement accuracy for H , 
operational data 

IX 

<7° 

°VV 

°VV 

1000 Km 

Efficient mode for taking operational data in- 
cluding wind speed along subtrack 

X 

0 

ff HH 

a 0 

°HH 

1000 Km 

Efficient mode for taking operational data in- 
cluding wind speed along subtrack 


Note: It may not be possible to get all-four polarization measurements made for all incidence angles at the instru- 

Cf 

mented test site during one satellite pass. 1 priority should be given to providing all incidence angles at 
each polarization, then to providing differential a 0 measurements at a given angle between polarizations. 


Figure 21.- Representative Mode Options for Scatterometer. 
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Figure 22.- Pictorial view of antenna beam illumination geometry. 



Location of beam edge after integration 



Figure 23.- Plan view of antenna illumination geometry (one side of satellite). 
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Figure 24.- Geometry of one resolution cell. 
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Note: V * Vertical polarization {transmit & receive) 

H - Horizontal polarization (transmit & receive) 


Figure 26.- Scatterometer mode 1 measurement sequence. 







Figure 27.- Transmit, receive, and processing timing diagram. 





Peak transmitter power, watts 


Altitude = 808 1cm 



Figure 28.- Transmitter power requirements for fan beam scatterometer. 
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Figure 31.** Margin available for various wind speeds. 


Normalized measurement error 
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Figure, 32.- Eadar scattering coefficient measurement error. 


Normalized measurement error 





-■* — — Sub-sateElite track - 

Figure 34.- Sample scan pa Hern el high resolution data cell illumination from forward and aft antenna -one polarization-node I. 
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Option n - Fixed increments of spacing (50 km) - variable incidence angles 
■ r = 50 km 


Beam illumination 
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Figure 35.- Geometry of surface illumination for one polarization - Mode n. 



Cell no. 15 14 13 12 11 10 9 8 7 6 a 4 3 1 



Cell no. 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 



Note - V = Vertical polarization 

H = Horizontal polarization 

Figure 36.- Mode II measurement sequence - two options 
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Figure 38.- Scatterometer block diagram. 









APPENDIX 1 

DERIVATION OF THE GEOMETRIC PARAMETER EQUATIONS 

By E. M. Bracalente 
Langley Research Center 

This appendix includes a description and derivation of the various equa- 
tions defining the geometric parameters involved in the scatterometer design. 
These parameters include angles and distances' defining the location of a data 
cell, dimensions of a data cell, and the doppler frequencies and bandwidths 
which form a doppler cell. 


Mode I 

Figures Al-1 and Al-2 illustrate the side view and plan view geometry 
associated with the fan beam illumination (also see figures 22 and 23 in the 
main text}. These figures should "be referred to for the derivation of the 
design equations. From figure Al-1, using the law of sines, we have that 

R 

sin 8. = :r— sin 6 (l ) 

l rv 


and 


R sin y 
c ” sin 0 


Using equation (l) and the relationship 

0. = 0 + y 

l 1 


the following equation is obtained, 


tan 0 - 


R sin y 
R^ — R cos y 


( 2 ) 


(3) 



The ground range, R , from the sub- satellite point is determined by the 

g 

relationship 


MEGEDING PAGE BLANK NOT FUAEEQ 


Al-l 


R = 


TT 


g 180 


■yE 


(5). 


Using these equations the location of each doppler cell with respect to the 
satellite can he obtained. 


The data cell locations along each fan beam is established by the proper 
choice of filter frequencies and bandwidths. To determine these filters 
certain requirements must be established. These requirements include resolu- 
tion (r), along-traek and cross-track spacing between data cells (S ) , swath 
width, and the approximate range of incidence angles. An additional require- 
ment is also imposed that two orthogonal data cell measurements must be made 
at the same surface of the earth. This is done using measurements from a 
forward and aft antenna which have azimuth angles that are 90 ° apart (see 
figure Al-2). 


The doppler filter frequencies, as shown a little later, are a function 

of the angle S, which in turn is a function of the earth angle y and 

ground range, R (see equations (4) and (5))< Since two orthogonal measure- 

g 

raents are required in each data cell, a value of y or R must be chosen so 

g 

that the aft antenna measurement overlaps the forward antenna measurement 

after the spacecraft moves the required ground distance (S^), (see figure Al-l). 

In this case, the value of & (data cell spacing along the antenna beam) 

c 

is chosen so that the cross-track spacing equals the along track spacing. In 
the following paragraph, the equations for the value of parameters S c , £ , 
and y that will satisfy these requirements are derived. The effects of 
earth rotation on the orthogonal measurement requirement and the doppler 
velocity are neglected here, but will be discussed in a separate section. 


The along-track and cross-track spacing S form legs of an equilateral 

S 

right triangle with the dimension 2& as the hypotenuse (for S ” 2r ). 

C 5 

The ground spacing R to the same instantaneous doppler cell on the forward 

g 

and aft antenna beam illumination form legs of an equilateral right triangle 
with the spacing S c as the hypotenuse (see figure Al-2). To a first 
approximation assuming the earth is a plane surface. 


Al-2 


( 6 ), 


& 

c 


2 


and 


c g loO 


(TV 


However, since these dimensions form a right spherical triangle on the surface 
of the earth, the following relationships from spherical trigonometry apply; 


cos(2& ) = cos^ S 
c s 


18 } 


and 


- 2 - 

cos S “ cos R 
c g 


(9), 


vftsre (25- ), S , S , and R are the earth center angles subtended by these 
esc g 

respective dimensions. Note that R^ is the same as the earth center angle, 
y , to the doppler cell. Re-arranging equations (8) ana (9) and substituting 
y in for R , the equations defining Z and S are obtained; 

££ C G 


, 18 0 S 2 

rt _ TlR "1 r St 

*0 “ 360 COB tcos 71 1 


do) 


_ ttR -1.2. 

S c = 180 COS (coB Tl 


( 11 ) 


The value of to the first data cell for an incidence angle of 

approximately 25° must be chosen so that proper overlap of the aft and forward 
data -.ell can occur. To do this S c must satisfy the equation 


= nS + S/2 
c Is s 


( 12 ) 


Al-3 


where n is an integer. The S / 2 factor in equation (12) correspond to 
J. s 

the distance traveled by the beam illumination during the scanning of the front 
antenna on each side of the S/C. The value of required for an incidence 

angle of 25° is approximately H6? km (808 l<m altitude), A value of n^ = 4 
will give an S value of 450 km for S = 100 km, and a value of n = 9 

C S X 

will give an S value of 475 km for S - 50 km. In both cases the incidence 
c s 

angle, 0^, falls approximately at the required 25 . The value of S c in 
equation (12) can be found (for a data cell at higher incidence angles) by 
increasing n^ by the factor (n-1), where n = the data pell number. Once 
the correct value of S c is found for cell number one, then proper overlap 
will occur (neglecting earth rotation effects) for all data cells when the 
spacing between doppler cells is placed at a distance £ . Note that only 
every other data cell is required to give equal along track and cross-track 
spacing if the spacing S is twice the required resolution. Substituting 

3 

equation (12) into equation (ll) and solving for Y, the equation for 
determining the value of Y^ (earth center angle to cell number l) is obtained 
as follows; 

Y ] _ = cos“ 1 [cos || S s (2n 1 + l)] 1/2 (13) 

For = 4 this equation reduces to 

810 S 1/2 

Y^ = cos _i [coB — R - ■ - -] (14) 


and for n - 9 to 

1710 S 1/2 

Y 1 = cos [cos — — ] (15) 

Equation (l4) is used (for Q08 km altitude) for S equal to approximately 
100 km (equation 3 in the main text), and (15) for S equal to approximately 
50 km. 

The next step in the design of the scatterometer system is the determina- 
tion of the doppler filter frequencies and bandwidths. Before this can be 


done the parameter L, which forms two sides of the Instantaneous Field of 

View {IFOV ) or doppler cell, must he determined. Figure Al-3 shows an enlarged 

view of a resolution data cell. As seen in this figure the IFOV is formed by 

the dimension L and the constant doppler lines (isodops) cutting across the 

heamwidth illumination determined by R 4>, where is the slant range from 

the satellite to the center of the IFOV and 0 is the antenna 3 db heamwidth. 

The parameter L is determined by the doppler frequency spread or bandwidth 

between the two constant doppler lines, f and f n . The subscripts 1 and 

u 

u stand for the lower and upper end of the IFOV. As pointed out in the text 

the idosops are hyperboloids on the surface of the earth (see figure 25 in 

the text). A line tangent to any point on this hyperbolic isodop will form an 

angle 6 with a line perpendicular with the line of S/C flight (see figure 

Al-3). The angle 6 increases as you proceed out along the beam illumination. 

Across the dimension R iji the isodop is assumed to be a straight line. The 

distance traveled by the IFOV during the measurement pei iod is equal to v t . 

The parameter t is related to the scan spacing S , the satellite ground 
P s 

speed v , and the number of antennas, (No^), and number of polarizations, 

(NOp), used in the measurements during one scan cycle. For the four antenna, 

two polarization case then. 


t 

P 


S 

s 


No a x N 0p x V g 



(16). 


The resolution data cell on the surface of the earth is determined by the area 
generated by the doppler cell resulting from satellite motion during the 
measurement period. The equivalent resolution cell is defined as the square 
whose dimension is the average of the dimensions of a rectangle within which 
the area generated by the moving doppler cell (IFOV) fits diagonally, as 
illustrated in figure Al-3. 

From figure Al-3 the equation defining the parameter L can be derived 
as follows. The equivalent resolution, r, is defined as 

r « (IT), 


Al-5 


■where 


a = c + J2. R <j> 


(18) 


and 


t = C + T t 

g P 


(19). 


Substituting equations (l8) and (19) into (IT ) and solving for c we have. 


c = 


2r - v t - /2Rd> 

■SB Si 


( 20 ). 


From figure Al~3 we see that parameter L = v^2d. Also from the figure we see 
that 


d = c - e 


( 21 ) 


where e = f sin 6 


and (from the law of sines) f = 


1/2 R <J) sin 135° 
c 

sin(l3f :C ’+ 6) 


R 0 

c T 


sin(l35°+ 6) 


( 22 ), 

(23). 


Combining equations (20), (21), (22), and (23), we obtain 


d = 


2r - v t - •J 2 . R $ R <j) sin 6 

SJ5... - 

2 sin(135°+ 5) 


'(au). 


from which we can write the equation defining L as 


v t /2 R (j) sin 6 

L = r R <j) - 

Si c 


sin(l35°+ <$) 


(25). 


All parameters in equation (25) are known or can be determined from 
previous equations except for the angle 6. The range R c used in this 
equation is the value calculated for the center of the IFOV (R^ varies 
slightly across the doppler cell from the f to f isodops). Also, since 

JJ 0 JJ 

* U 
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the angle 6 varies slightly across the doppler cell, the value of 6 at the 
center of the cell is used. Before the equation defining 6 is derived, the 
equation defining the doppler frequencies and hand-widths must he derived. 

Figure Al-h is a pictorial view of the antenna heam illumination along 
with appropriate velocity vectors. The earth velocity vector orientations are 
shown for the equator crossing for both the ascending and descending direction. 
The earth velocity vector will he neglected for this first derivative of the 
doppler frequency equation. A later section will cover the effects of earth 
rotation on the doppler frequency. By inspection we see that 


and 


v. = v cos d) * 
Is T 


( 26 ), 


v.^ sin 0 


(27), 


where v = the satellite velocity 
& 

v 1 = the component of satellite velocity in the plane of the fan heam 
and v^ = the component of doppler velocity or range rate along the 
range line R c . 

The two-way doppler frequency is found from the well known equation 




2 v,f. 
d t 


( 28 ), 


whe e - the transmitting frequency, 
and c = the speed of light. 

Positive doppler frequencies occur on the two antennas pointed forward of the 
satellite and negative doppler frequencies on the antennas pointed aft. 

Comhining equations (26} , (27), and (28) the following equation results; 

2f 

fp " ± Vs cos sin 0 ( 29 ). 


Since the antennas are placed at b 5° to the satellite velocity vector, 
equation (29) reduces to 
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£ x + 

D — c 


v sin 8 
s 


( 30 ). 


To determine the doppler bandwidth, the doppler frequencies f and f_ 

4J ^ D 

for the XFOV must he found. Shis is done hy first determining R and R u 

S ^ Sy, 

from the equations 


R = R - L/2 


g- 


S, 


(31) 


'n 


and 


R = R + L/2 
s u S n 


(32) » 


where R and R are the surface distances from the subsatellite point to 

the upper and lower ends of the XFOV and R is the surface distance to the 

Sv 


center of the XFOV for cell number n. R a is determined from the equation 


g. 


n 


R « R + (n ~ l)Jl 


g n S 1 


(33) 


where R is determined hy equations (13) and (5)i & c I s determined by (10), 
®1 

and n is the doppler cell number along the beam illumination. Using the 

values of R and R , the angles 0 and 0 ft can be determined using 
61 g ^ u & 

equations (4) and (5). The dopper bandwidth B c is then found from 


B = f 
c D 


- f . 


u 


Or 


(34): 


Using equation (30) this becomes 


& f , v 

B c = “ [Bin 0 U - sin 0^3 


( 35 ). 


Having derived the doppler frequency equation we can now derive the 
equations needed to determine the value of 6. This can he done hy solving 
for the angle ifj (OCA) in triangle AOC of figure Al-5s since by inspection 
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6 = ^5° - ^ (edge of "beam illumination assumed parallel to beam centerline). 
It is assumed in this derivation that the lines OB and OC are coincident 
(this is not quite the case* but the deviation from coincidence is s m all and 
can be neglected.) Using an equation of the same form as equation (U) and 
triangle AOCjthe equation for t|J is written as 


R 


\p'~ tan 


-1 


6, 


sin 4) 


n 


g 


x 


6. 


COS (f) 


n 


S 


(361. 


To solve for the angle iJj the dimensions R , R and <J> must he 

determined. R (ground range for cell number n) is determined from 
§ n 

equation (33). 0 is determined from the right triangle OAB by the 

S 

relationship 

R <f> 

*g = tan ' 1[ ri _] 


g 


(37), 


n 


where R^ is found for cell n from equation (2), and (j> is the antenna 
3db beamwSdth in radians . To determine R^ , the doppler frequency f ^ 
for the center of cell n is first calculated . Since the doppler frequency 
is the same anywhere along line DC, the angle B is found using equation ( 30 ) 


, c fn 

„ . -I D n 

0 = sm r 

x 2 v f . cos 0 


(38) 


s t- 


X 


where (j> = tj)' + (see figure A1-.5). 


r x 


The dimension R is calculated using equations (l), (3) and (5) to obtain 


g. 


x 


R = R[(sin‘ 
g„ 


-1 A 


R. sin 0 


R 


x - V ife] 


(39). 


how solving for the angle ^ we can obtain the angle 6 from 


6 = 45 - ifi 


(fco). 
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Mode II 


Equations (l) through (5) are applicable for the mode II design. Figure 

35 in the main text illustrates the geometry of the surface illumination for 

mode II. Also included in this figure is an enlarged view of the resolution 

data cell. The along track dimension of the resolution data cell is deter- 

mined partly "by the doppler bandwidth, which establishes the parameter L, and 

partly by the smear distance generated during the measurement period. The 

cross track dimension is formed by R c (j>. Since the antenna beam for mode II 

is directed along the antenna flight line, the resolution data cell generated 

is long and narrow. The resolution in this case is defined by the sum of L 

and the distance v t traveled during the measurement period, t . Any 

g p P 

combination of doppler cell length (L) and measurement period (t ) can be chosen 

k 

to form the specified resolution length. As shown in appendix 2 it is 

desirable to maximize the product B where B is the doppler bandwidth, 

c SI c 

is the net integration time for signal plus noise during each 


and 


SB 


measurement period. is a direct function of the parameter L, and T 

is a function of t . For a fixed resolution, r, there is optimum choice of 

P 

L and t which will maximize S T_„. The following derivation shows that 
p c SB 

this occurs when tv = L. 

P g 

To determine the optimum choice of L and t^., an equation for 
is written in terms of L and r. The derivative of this expression with 
respect to L is taken and set equal to zero. The resulting equation is then 
solved for L. This is shown as follows. By definition 


or 


r - L + t v 
P S 


r - L 


p v 


(hi), 

(42); 


g 


Also 


T = — T 

SB T L P 


( 43 ), 


where T = the pulse repetition period of the transmitted rf pulse, and 
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Tp s= the pulse width of the transmitted rf pulse. Substituting (42) into 
(43) we have that 


(r - L}t p 

T SU = v“T 
g 


(44). 


Since at any incidence angle the bandwidth B c is linearly related to the 
dimension L, for small changes in L, the following proportionality can be 
written: 


B 

< 

B 


c L 


(45), 


or 


B - B f- 
c c x L x 


( 1 * 6 ) 


The equation for B t„. t is now written as 

c oil 


B c T SN 


B t (r - L )h 
C 1 P 
L i T 

1 g 


(47). 


Combining all constants and setting them equal to k, equation C4?3 is written as 


B T_,. t = k L r - k If 
c SN 


(48), 


To maximise R we take the derivative of equation (48) and set it equal 

o SN 

to 0; thus. 


a < E cV 

dL 


= 0 = k r - 2 k L 


(49) 


from which 


r - 2 L 


(50). 


Substituting this value of r into equation (41) gives the result 
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( 51 ) 


Tims we see to optimize the mode II design for a fixed antenna gain and 
transmitter power, the measurement period and bandwidth must be chosen so 

the IFOV length h , equals the integration smear distance, v t . 

6 P 

Since the antenna beam is directed along the satellite direction of 
flight, the angle cj>' in equation (29) is zero* Therefore the doppler 
equation used in mode II calculation becomes 

2 f 

fp = + v sin 0 (52) 

— c s 

Similarly the equation for calculating the doppler bandwidth for <j)* equal to 
zero degrees, is from equation (35) 

2 f v 

■q = - j Bin q _ sin e ] (53) 

c c u x, 

Earth Rotational Effects 

In the final system design, the earth rotational surface velocity must 
be included along with the satellite velocity in determining the doppler 
filter frequencies and bandwidths. Figure Al-*6 shows a plan view of the 
satellite track on the surface of the earth along with pictorial views of the 
beam illuminations at different parts of the orbit. As seen from this figure 
the orientation of the satellite velocity vector and the antenna beams in 
respect to the earth velocity vector does not remain the same throughout the 
orbit. At the peak north and south latitudes, the satellite velocity vector 
is l8o° away from the earth velocity for the retrograde orbit. Throughout 
the remaining portions of the orbit the satellite velocity vector and antenna 
beams rotate + J2° from this position for a 108° inclined orbit. The 72° 
orientation occurs at the equator crossings, as shown in figure Al-6. 

In the previous derivation of the doppler frequency equation (equation 
30) the earth surface velocity was not included. In figure Al-4 the 
orientation of the earth velocity vector is shown, at the equator, for both 
the ascending (south to north travel) and descending (north to south travel) 
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directions of the orbit. Note that the earth Telocity is shown in the plane 
of the satellite translated from the doppler cell. Mote also that since the 
projected antenna "beams intercept the earth away from the sub-satellite point, 
the satellite and any particular doppler cell are not at the exact same 
latitude and longitude. This fact will be neglected in the following 
derivation of the doppler velocity equation that includes earth velocity. 

The effects of the latitude and longitude differences will be discussed later. 

From equation (26) we have that v. ” v cos <f? 1 . The earth velocity 

.L *5 

component (v^) in the plane of the antenna beam must be combined with v^ to 
determine the total doppler velocity v^. From figure Al-4 we see that for 
beams 2 and 4, 


and 


cos ^ (ascending) 

(54) 

cos (descending) 

(55) 


4 


where v = earth rotational surface velocity, 

6 

and and = the angle between the earth velocity and the plane of 

the antenna beam for the aseending and descending directions respectively. 
The earth surface velocity is a maximum at the equator and decreases with 
increasing latitude, with a value of zero at the poles. The earth surface 
velocity at any latitude can be shown to be 


v = V cos <J) T (56) 

e e L 

where d> T is the earth latitude angle, and V is the earth surface Velocity 
(.465 km/s ) at the equator. The doppler velocity v^ to each doppler cell 
is found from the vector sum of v^ and v^ multiplied by the sin 0. 
Combining equations (26), (54)* (55) and ( 56 ) we have for antenna beams 2 
and 4, 
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and 


v, = +[v cos * + V cos <j> cos tJj. ] sin 9 (ascending) 

CL S 6 ij a 


( 57 ) 


v. - +[v cos <j)’ - V cos <L cos Tk,] sin 6 (descending) (58), 

ct “ s e jj u 

where the positiYe value corresponds to the forward antenna ("beam 4) and the 
negative value corresponds to the aft antenna (heam 2), Similarly equations 
(57) and (58) also apply to beams 1 and 3 except that the signs in the brackets 
are reversed. Referring u ek to figure Al-4 we can see for = 4l? that 
angles 


and 


*A = “i " 45 ° 


ift D = 135°- to. 


(59) 

(60) 


for beams 2 and 4. Similarly for beams 1 and 3, 


and 


* 135 - “i 


= 0 ). - 45 


(61) 

(62). 


The term in these equations, is the angle between the earth’s surface 

velocity and the negative satellite velocity. Equations (59) and (6l) 

apply when v is iu the 4tn quadrant (ascending) and equations (60) and 
e 

(62) apply when V g is in the 3rd quadrant (descending), where quadrants 
are numbered counter-clockwise from beam number 1. As mentioned before^the 
angle varies + 72° throughout the orbit. At any point in the orbitj the 

angle Oh can he shown to be a function of the orbit inclination angle CL 
and the latitude position (angle cj) ) of the satellite. The equation 
defining the angle Oh , which can be derived by vector analysis , is given as 


, cos a. 

—l r 3. n 

03. « cas [- - ■■ 1 " 7 — ] 
1 cos <p T 


• (63). 


Combining equations (59) or (6l) with (57) and equation (60) or (62) with (58). 
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I A 

along with the proper sign choice inside the bracket, and for <£' = 45 , the 
following two equations result which are used to calculate the doppler - 
velocities on each antenna beam. 

v^ = +[^| v s + v e cos cos (ok - 1*5°)] sin 0 

j and 

I J? 

‘ v- = +[ — •x v - V cos <f> T cos(l35°- w. )]sin G 

j CL — * d. S 6 ii 1 

[ 

I 

j where un is found from equation (63). In the ascending portion of the orbit, 

| equation (64) is used for beams 2 and 4, and equation (65) for beams 1 and 3- 

j The equations are interchanged for the descending portion of the orbit. 

Substituting equations (64) or ( 65 ) into equation (28) results in the equation 

! 

I for doppler frequency. The equation for the bandwidth, B p , is obtained by 

] substituting in the resulting expression for f^ into equation (34). 

Although a secondary effect in the overall doppler velocity equation, 
j the difference in latitude and longitude location between each doppler cell 

j 

1 and the satellite will have to be considered when determining the actual doppler 

i 

i frequency for each data cell (or location of data cell, given doppler 

I frequency). Equations (64), and ( 65 ) will have to be modified to account for 

: this effect. In these two equations , the latitude angle $ of the 

particular doppler cell must be used, rather than the latitude angle of the 
satellite. Also, the second term in the brackets of equations (64) and ( 65 ) 
must be multiplied by a cosine term, which is a function of the longitude 

\ difference between the satellite and the particular doppler cell. This 

| 

I correction can be seen by a study of figures Al-4 and Al-6, As seen in 

| - figure AIj- 4, the earth velocity translated from the sub-satellite point is 

| in the same plane — which is parallel to the surface of the earth — as 

i 

1 - the satellite velocity vector. However, when translated from any particular 

1 

doppler cell, the earth velocity vector will be rotated out of this plane, 
due to the earth curvature. This rotation occurs because the earth surface 
velocity vector is always perpendicular to a plane containing the surface 
point and the axis of rotation. It can be seen that the angle of rotation 
for the scatterometer case is the difference between longitudes of the 


(64) , 

(65) , 
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satellite and the doppler cell. The second term in the brackets of equations 
(64) and ( 65 ) must be multiplied by the cosine of this angle. Hence, 
equations (64) and { 65 ) can no w be written as 

V d = ±[\£T /2 v s + V g cos 4> cos cos(tii - 45) 3 sin 0 (Mode I) (66), 

V d " ±J^ 2 ./ 2 v g - V g cos (j) cos cos(l35°” tin )]sin 8 (Mode I) ( 67 ), 

where is the longitude difference between the sub«satellite point and 

a doppler cell n, and <L is the latitude of donpler cell n. 

Ln 

By similar analysis and by referring to figure Al-4, the equation for 
the total doppler velocity on doppler cells in Mode II can be written as, 

V. - + [v + V cos 0 cos AQ cos to.] sin 0 (Mode II ) (68), 


The plus sign corresponds to the antenna directed forward of the satellite and 
the minus sign for the antenna directed aft of the satellite. 

The term in the brackets of equations (66) and (67) is the 
component of doppler velocity in the plane of the antenna beam parallel 
to the surface of the earth. This is the velocity v^ (ascending) shown in 
figure Al-4. As seen from equations (66) and (67) s the absolute value of 
the velocity component v^ will vary in opposite directions on each antenna 
beam throughout the orbit. The variation of this velocity component for 
each antenna beam is shown in figure Al-7 (effect of latitude and longitude 
difference neglected in this figure). As seen from this figure the mean 
of the absolute value of v^ occurs at the northern and southern latitudes 
of 72°. At this point in the orbit the value of v^ is the same on each beam 
(neglecting the differences that occur due to the difference in the latitude 
and longitude location of each doppler cell). The peak excursion from the 
mean value occurs at the equator. The doppler velocity component v^ will 
vary approximately +6% from the mean value (or total of 12%) throughout the 
orbit (see figure Al-7). The result of this variation will he to shift the 
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location of the data cells in and out along each "beam, since only one hank of 
fixed doppler filters will he provided in the spacecraft. The mean value of 
velocity v^ will he used to calculate the doppler filter frequencies and 
handwidths. An additional variation of less then 1% on the outer doppler 
cells and less than .1% on the inner cells mil occur, due to the effect of 
the latitude and longitude difference between the satellite and the doppler 
cells . 

The variation of the velocity v^ could be made small (within the 
variation caused by the latitude and longitude differences} by rotating the 
spacecraft or the antennas hy a small amount around the yaw axis, as a 
function of spacecraft latitude. The effect of this yawing is to compensate 
for the variable earth velocity vector by varying the component of satellite 
velocity -tfhich contributes to the total doppler velocity. The amount of 
spacecraft rotation required at each latitude can be found by re— writting the 
terms in the brackets of equations (64) and (65) as 

v_ = + [v eos(45 + A) + V cos cos(u). - 45°^ A)] (69) 

3 — S 6 L 1 

and 

v = + [v_ cos (45 + A) - V cos <j> T cos(l35 + A - w. )] (TO) 

Setting them equal to each other, and solving for the term A - where A is 
the amount the S/C or antennas are rotated - in terms of <f>^, tit , and the 
satellite and earth velocities , the following equation for A results} 


cos cjn sin gj. 

tan A = + — -XZ — — t~~ — — ~~~ (71 ). 

— V /V + COS tp_ COS tlJ. 

s e L i 

The positive sign is used in the descending portion of the orbit and the 
negative sign in the ascending portion. Solutions for A in equation (71 )» 
at various latitudes throughout the orbit , show that a maximum rotation require- 
ment of approximately +3.3° occurs at the equator. The value of rotation 
required decreases with increasing latitude. Zero rotation is required at the 
northern and southern extremities of the orbit. Substituting the value of 
A into equations (69) and (70 ) for each latitude shows that less than .2% 
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variation in the velocity v^ occurs throughout the complete orbit. 

As pointed out in the main text, the orthogonal measurement capability of 
the scatterometer system is greatly affected by the earth rotation. The ortho- 
gonal measurements on a particular surface of the earth are provided by the 
forward and aft antennas. However, between these measurements, the aft beam 
illumination must move the distance S c , the instantaneous surface distance 
between the forward and aft doppler cell locations. For the baseline design 
this distance is 1*50 Ion for cell number 1 and 1850 km for cell number 15 (see 
figure 29). For the ground speed of 6.6l km/s, an elapse time of 68 seconds 
will occur between the orthogonal measurements of data cells number 1 and 280 
seconds for data cells number 15. At the equator, during this elapse time 
period, the surface area measured by the front antenna for cells number 1 and 
15 will move 31 and 128 km respectively due to earth rotation. (The 
remaining data cell measurements will have moved an amount between these two 
extremes). As the spacecraft travels north or south from the - equator the 
amount of surface movement decreases since the earth surface velocity is 
decreasing. Figure Al-8 illustrates this effect for various data cells at 
different latitudes. In this figure the solid and dotted squares represent 
the location of the aft and forward antenna data cells, respectively, for a 
given orthogonal measurement pair. The arrows drawn between the aft and for- 
ward measurement cells represent the direction and magnitude of movement 
caused by earth rotation. 


The surface distance, S^, separating the orthogonal measurements for any 
set of data cells is found from the equation 


S 


m 



cos <j> 


Ln 


(72) 


The spacing, S , for each data cell is found from equation (ll). The angular 

direction of movement of the front antenna data cell, measured from the 

m 

satellite subtrack, is approximately 


* ffi - i (iso 0 - «i > ± an* 


(73) 
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The positive signs are used during the descending portion of the orbit, with 
counter-clockwise rotation positive, looking down on the earth. Pages U5 and 
1+6 of the main text discuss the effect of nonoverlap of the orthogonal mea- 
surements, and recommends a possible solution to this problem. 

Before the doppler velocities from equations (66) and (67 ) can "be calcu- 
lated! the equations for the latitude of the doppler cell (tj)^) and the 

longitude difference (Aft ) between the satellite and doppler cell must be 

hn 

determined. The equations for these parameters are derived using figures Al-9 
and Al-10. The equations for the latitude and longitude differences at the 
equator are shown in figure Al-9* They are determined from spherical trigono- 
metric relationships for right spherical triangles. For small Y n (<15°)s the 
equations in figure Al-9 can he simplified to the following equations: 

Beams 1 and 3 


Afi = 
Ln 


sin(w. 

i 

+ 45°)] 

m) 

A hn = 


cos (w. 
X 

+ 45°)| 

(75) 


Beams 2 and 4 


A n = 

Ln 

It 

1 'n 

cos(oj^ + 

45°) | 

(76) 

ll 

a 

-o- 

<1 

1 T 
' 'n 

sin(ti3. + 
i 

45° )| 

(77) 


where Aft^ and Ai})^ are the longitude and latitude differenceSj 

respectively. The minus sign inside the brackets is used for the ascending 

portion of the orhit and the plus sign for the descending portion. These 

equations are exact only at the equator. However, for small y n angles the 

latitude of a doppler cell (<j) ) can be approximated by the equation 

Ln 




+ Atj). 


Ln 


( 78 ) 


where A6_ is found using equations (75) and (77)* Also it can he shown 
Ln 
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that can be approximated at any longitude by dividing equations (7*0 and 

(76) by the cob <j) . The longitude angle, of a doppler cell is found from 


SL = £2_ + A£L 
Tin L — hn 


(79). 


Note, that only the magnitude of A<{) and Afi^j found from their defining 
equations, should be used in equations (73), (78), and (79)- The choice of 
sign used in equations (78) and (79) is a function of the relative location of 
the satellite and the doppler cell with respect to the equator and the zero 
meridian, respectively. 

For the larger incidence angles (outer doppler cells) at the higher 
latitudes the exact equations defining (j> and A& T ^ may be required to 
determine' the doppler cell locations with sufficient accuracy. The exact 
equations defining 4>-^ n and A£J^ n are found from figure Al-10. Using the 
shaded spherical triangle in figure Al-10 and the law of cosines for sides 


" sin ^tsin d) T cos 
Ln , T L 


Y + cos d>_ sin 
n L 


y cos a] 
J n 


( 80 ), 


and using the law of sines 


sin a sin y 

ASL = sin -1 [ t ~] 

Ln cos $ Ln 


( 81 ). 


where, in the northern hemisphere, a is the angle measured from north to the 

antenna beam center (measured from south in the southern hemisphere). 

The angle a is a function of the angle a), (angle between the satellite 

1 

and earth velocity vectors). The equation defining a is found from one 
of the following equations; 


a = au + 45° (82) 

a = 135°+ w i (83) 

a = 225°- a> ± (84) 
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where the equation choice is based on the hemisphere location of the satellite 
and the direction of the orbit (descending or ascending). TabDe Al~2 lists 
the equation choice for each antenna beam for all portions of the orbit . 

Table Al-1 summarizes the equations which have been developed for use in 
designing the scatterometer . As noted before, the doppler frequency at a 
specific incidence angle varies throughout the orbit. The mean doppler fre- 
quency occurs at the maximum northern and southern latitudes of the satellite 
orbit, as seen in figure Al-7. However, since the location of an antenna 
beam changes latitude position with respect to the satellite at the maximum 
satellite latitudes (see figure Al- 6 ), the mean doppler frequency at any 
particular incidence angle will be slightly different at these two latitudes 
extremes. Therefore, for each incidence angle, the doppler frequencies and 
bandwidth^ are determined by the average of the mean value calculated at the 
maximum northern and southern latitudes. This procedure is used in section II 
of table Al-1 for determining the doppler filter frequencies. 

Table Al-2 lists the pertinent equations used for determining the doppler 
cell location and the parameter L (for use in eq. (5) of the main text ) , 
given the doppler filter frequencies and bandwidths. Equation ( 6 ) of table 
Al-2, which determines the angle is a transcendental equation which must 

be evaluated by a digital computer in an iterative process. Three additional 
equations s imi lar to equation ( 6 ) can be written for beams 1 and 3 in the 
ascending orbit and for beams 2 and 4 in the ascending and descending orbit. 
The exact equations for the latitude angle, <J) T , and the longitude difference 
angle, Aft^, are listed in table Al-2. The approximate equation for these 
parameters ( 78 , 74, 76 ) could be used for the smaller incidence angles at the 
lower latitudes, thus simplifying to some extent equation ( 6 ), of table Al-2. 
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STABLE Al-1 SOATIEROMETER DESIGN EQUATIONS 


I. Geometric Equation of Doppler Cell Location 
-1 R A 

1. 0. = sin sin 0 

X Jtv 


2 r = R sin 
c ” sin 0 


,, o _ ^ -If R sxn Y i 

3. 0 = tan L" = — — ~i 

R^ - R cos Y 


k. R - Y R , Y in radians 

g 

, 180 S 2 

r- n ?VR “1 r Si 

5 - £ c = 360 005 [cos “IS ' 


6. E = B + (n - l)t 
s n E 1 c 

T- Y x = cos~ 1 [ cos s s + l).] 1/2 


8. n = S /S e - 1/2 , n. is nearest whole number to 

X C S X 

satisfy this equation 


Equation 
number 
in text _ 


( 1 ) 

( 2 ) 

(M 

(5) 

do) 

(33) 

(13) 

( 12 ) 


9 * S c “ lfo eos ~ 1 ( cos2 Y ) , for e t = 25° 


( 11 ) 


II. Equations for Determining Doppler Filter Frequencies 


1* f D 

n 


f D (N) + f (S) 
n n . 


where f„ (JJ) = Doppler frequency for cell n at maximum northern lat. 
■n 


and f^ (S) = Doppler frequency for cell n at maximum southern lat. 

"n 

1 / 2 . f 

2. - ± — n — £ v + V cos cos AJ2 ]sin 0, Mode I «f>- — max. (28), (66) 

n C S e Ln Ln L or (6?) 
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TABLE Al-1 (continued) 


3* - sin ^[sin ^ cos y^ +_ — - cos ({j^ sin y^] , (f)^ = maximum (8o) , (82) , (83) 

+ maximum northern latitude 
- maximum southern latitude 

& sin Y n 

lu A SL --sin [ _ ■ — r— 1 , (j) T = maximum 

Ln L 2 cos <p, ’ ^L 


J2 


<ai) 


Ln 


2 f. 


5. f, - + ^ [v + V cos cL cos A£2_ cos to. 3 sin 0 , Mode II (68) 

Lf ~~ v# S 0 Lrll i-lll 1 

n 

See Table Al-II for equations defining, «(> T , AQ^, and on for equation (5). 


III. Equations for Calculating Parameter L 


^r) u xu x v 
F OA OP g 


8 v » 

s 


for 


H antennas 
2 polarizations 


(16) 


R if) 

2. 4 = tan' 1 [jf-] 


g. 


n 


_ - . -1 
3* sin 


C f. 


n 


jr ■ 
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(28), (27) 

( 66 } 


Note: Use the average values of tj}^ and A£L ^ calculated for the maximum 


Ln 


northern and southern latitudes, use equations II- 3 and. II-U of this table. 
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TABLE Al-1 (concluded). 


6 . 6 = 45 - ^ 

V t v'cf R (jj sin 6 

T. L = ^r--^-R^ . £■-- v f 

^ c v sin (135 + o) 
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IV. Equations for Calculating Doppler Bandwidth, B 
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s i S n 

2. R = R + L/2 

S u S n 
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g n g l 
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ft 
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R 
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(31) 

(32) 

(33) 


00. (5), 
(31) 


00. (5), 

(32) 


6. B = f_. - f _ = * — ~ - [v + V cos (j> cos AJL ][sin 0 - sin 0„) 

c D u C s e Ln Ln u V 


4> l = max. 


( 28 ), ( 66 ) 


Note: Use the average values of 6_ and AS2 r calculated for the maximum 

Ln Ln 

northern and southern latitudes, use equations XI- 3 and II-4 of this table. 


TABLE Al-2 EQUATIONS FOR LOCATING DOPPLER CELLS 


Equation 
number 
in text 


2 f r 

1. I - - ■ — [--| v + V cos 4>_ cos A£L cos(to. - ^5°)3sin 0 (66), (28) 

1 D ‘ C 2 s e T Ln x n 

n 

Ascending - Beams 2 and 4 
Descending - Beams 1 and 3 


2 f >- 

If I ^ — [-^ v - V cos 4>r cos ASI. cos{ 135 °- to,)]sin 0 

1 D * C 2 s e Ln it in 

n 

Ascending - Beams 1 and 3 
Descending - Beams 2 and 4 


3- 

k. 

5. 


-1 


cos a. 


cl = cos [ j“] 


"COS ^ 

“1 T 

d» T = sin [sin tJ> T cos y + cos <Jj sin y cos aj 
Ln L n L n 


_ sm a sin y 

* -lr Ht 

ASL = sm v J 

Ln cos <}) 


r Ln 


(67), (28) 


(63) 

( 80 ) 

(81) 


a Determined as Follows 


Northern Hemisphere 

^ , ,,-0 + Desc. 

Eeam 1 a = to. +. 45 - 

„ „ , . _ o + Asc . 

Beam 2 a = <4. + 45 - _ Desc . 

Beam 3 a = 135°“ CL - DESC 

a = 225°- to. - ASC, to >45° 

1 -l* 

a = 135°+ a - ASC, cl < 45° 


Southern Hemisphere 

Beam 1 ct “ 225°~ CL - ASC, an > 45° 

a = 135°+ to. - ASC, to i < 45° 

a = 135°- cl - DESC 

Beam 2 a = 135°“ \ ~ ASC 

a = 225°- to ± - DESC, ol > 45° 

a = 135°+ cl - DESC, cl < 45° 


Al-25 


TABLE Al-2 (concluded) 


Beam 4 a *= 135 to, - ASC. 

x 

a = 225 - to. - DESC, w. > 45 
x x 


Beam 3 a - ta. + kf * ® esc * 

1 “ aSC • 


a = 135°+ 0 ). - DESC, W. <45° Beam 4 a = ax + 1*5 ° + f 80 ’ 
xx i “ jjssc • 


Mode II a - 90° - ax - ASC 

1 ~~ * • x 

a = 90° + a) - DESC 


Mode II a = g0° + ax - ASC 
a = 90° - ax - DESC 


D R sin \ JS 

TfT = Bin[tan _ ( - _ ■ >K-g v + ? e cost Bln" (sis * cos Y n 

t A ' n 


+ cos cfa T sin Y cos a)] x 
T L n 


cos [sin -1 (- 


sin a sin y 

'n 


cos[sxn ^(sin tf> cos Y + cos tp sin x, cos a)] 

• L n n 

cos (ax - 45)) 

—1 

Equation 6 is formed from equations 1, 4 and 5 of this table, and 0 - tan 

R sin Y 

[- B ~' d - equation 1-3 of table Al-1. Equation 6 is used for beams 1 

R. — R cos Y 
xxr n 

and 3 in the descending orbit. Similar equations can be written for beams 1 
and 3 in the ascending orbits and for beams 2 and 4 in the ascending and de- 
scending orbit, ox is found from equation 3 of this table, a is chosen from 
the group of equations on the previous page of this table. 

The parameters 0., R and R can be found from equations 1-1, 1-2, 

X C ■££ 

n & n 

and 1-4 respectively of table Al-1. 


T- 


8. L = 2y^ R - .2 Rg , y^ in radians (Equation 32) 


9- Y = sin _1 (R./R sin 8 ) - 3 (equations 1 and 3) 


u u 


10. 0 is found from equations 1 or 2 of this table where |f_ ] is 

u D H 

replaced with f_ 
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Figure Al-2.- Plan view geometry. 





For beams 2 and 4 
v 3 = + (Vj + v 2 )(ASC) 

v 3 = + (vj - v 0 )(DESC) 

\p A - u). - (b 1 = u>. - 45° 

l ^ l 

if/jy = 180 - 4>' - w. = 135 - w. 
For beams 1 and 3 
v 3 = + (vj - v 2 )(ASC) 
v 3 = ± ( v i + v 2 )(DESC) 
j^ a = 135 - w. 

= W i ‘ 45 

Note: To determine if earth 
velocity adds to or subtracts 
from the satellite velocity, 
place earth velocity vector 
at the Doppler cell 
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Figure Al-4.- Pictorial view of antenna beam illumination showing satellite and 
earth velocity relationship at the equator. 




Figure Al-6.- Plan view of satellite track showing pictorial views of the 
beam illuminations. 



Magnitude of earth and satellite velocity component, v 3 , km/s 




Figure A1 -7.- Variation of velocity component, v^, throughout the orbit. 


ORIGINAL PAGE IB 
OF POOR QUALTTS 




Figure Al-8.- Relative locations of forward and aft data cells (beams 1 and 2) when 
orthogonal measurement is made - for three satellite latitude 
locations, ascending portion of orbit. 



( 1 ) 



Figure Al-9.- Calculation of latitude and longitude difference. 



used to find latitude (c^ n ) 


APPENDIX 2 


DERIVATION OF LINK AND MEASUREMENT ACCURACY EQUATIONS 

By E, M. Bracalente 

Langley Research Center 

Link Equations 

The equation (equation (5) in the main text} which determines the value 
of the radar return signal power can he derived as follows. The power density, 
P^, existing on the earth's surface at a doppler cell is 


F 


a 


p t ° 

i -2 

4 ir R 

c 


(1) 


where P^ is the transmitter power, G is the antenna gain to the center of 
the doppler cell and R c is the slant range from the satellite to the doppler 
cell. Since the transmitted pulse length is relatively long, the antenna heam 
is completely filled with the signal. Therefore, the total power radiated hack 
toward the spacecraft from a doppler cell is 


P. K o 
d e 


P t ° A c g ° 

k TT R 2 
c 


(2) 


where A q is the doppler cell area on the surface of the earth and 0° is the 
normalized radar scattering coefficient. The power density at the spacecraft 

O 

is equal to equation (2) divided hy 4 ttR c . The total power at the receiver 
input is then equal to the product of this power density, the ant enna 's 
effective collecting aperture, and any additional system losses, i.e.. 


P R - 


_VL 

It It r! 


A 0° #4 

rt ' fl nr' 


4 7T , 2 

4 TT R 

c 


(3) 


where 

GX 2 

- antenna effective collecting aperture, 

L = miscellaneous system losses (ie. , two way atmospheric, receiver, and 
transmitter losses), 

A^ ™ R^L “ doppler cell area. 
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G 

G = G o {^“) = antenna gain to the center of the doppler cell, 

Utt e 

G « = peak antenna gain, 

O <j)g 

d> = 3dB heamwidth across beam, 

3 = 3dB beamwidth. along fan beam, 

L = Doppler cell length, 
e = overall antenna efficiency. 

Substituting these relationships into equation (3) and rearranging, the follow- 
ing equation for the received power is obtained. 


P R " 


F (G/G ) 2 0° A 2 e 2 L L 
t o s 

Utt R 3 0 6 2 
c 


(M- 


The received signal-to-noise ratio from a doppler cell is. 


S/K = 


'R 


K B 
s c 


( 5 )» 


where T s is the system noise temperature and B c is the I.F, bandwidth which 
equals the doppler spread bandwidth across one doppler cell. The system noise 
temperature, T g , can be referred to either the input to the rf preamplifier or 
at the antenna terminals. Either place is correct, but if T g is referred to 
the preamplifier input, the loss factor L„ must be included in the link 

L\ 

equation. If T g is referred to the antenna terminals, the factor is not 

included in the link equation. Figure A2-1 shows the derivation of the equa- 
tions for both of these cases. 

Derivation of the parameter =*> 

As shown in figure A2-2, a gated integrator is used to average the output 
of a square lav detector over a number of return radar pulses. This measure- 
ment P is an estimate of the mean return signal plus noise power. Next, a 
measurement of the antenna and receiver mean noise power P^ is made by 
integrating the square law detector output in the absence of the radar return 
pulse. P^, is obtained from these quantities by 

k\ 


A2-2 


( 6 ), 


P = P - P 
r R SN If 


If 


M 


P gH and are assumed to "be independent, then the standard deviation 


of P, 


is 


°M = *-°SK + 


1/2 


( 7 > 


where 


^ Var P SN^ 


the variance of the measurement of signal plus noise, 
and 


SEf 

2 , /n . 

= (Var P^) = the variance of the measurement of noise only. 

2 

Referring to figure A2-2, the equations for the variances and a. 


33 


'if 


can "be determined as follows. Davenport and Root shows that a Gaussian 


signal uniformly distributed over a bandwidth B c , has a triangular spectrum 
at the output of a square law detector and low pass filter. As illustrated in 


figure A2-2, the peak value of this distribution is where 


is the 
u 

one sided power spectral density of the signal in the bandwidth B q . The 
distribution decreases linearly to zero at a frequency of _+B , from this peak. 
The mean value of the signal is equal to . 

The radar return signal is assumed to be a Gaussian random process uniform- 
ly distributed over the bandwidth B c . The power density for this signal is 
P R /B c . The power density of the antenna and receiver noise, which is also 


uniformly distributed and Gaussian, is 
For signal plus noise, then, 


KT g , 


where K is Boltzman's constant. 


and for noise only, 


P D = ET S + P k /B o , 


P D “ ^ 


( 8 ) 


( 9 ). 


The spectrum at the output of the gated integrator is shown in figure A2-2. 
It has a (sinx/x) ^distribution. The variance, a P , at the output of the 
integrator is equal bo the integral of the (sinx/x ^function times the input 
spectrum, i.e. , 
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B 

o = 2 f Q ° S N (f)H(f)df 


( 10 ) 


where ( 1* ) is the spectrum input to the integrator. 


S„(f) = P^(B rt - f), (see figure A2-2! 

JN DC 


( 11 ) 


and H(f) is the power transfer function of the gated or running mean integrator. 


Htf) = sid (f . n 

(«r 


( 12 ), 


where T is the total integration time. Substituting equations (8), (ll), 

2 

and (12) into (10), the variance a gfI for signal plus noise is 


B 


• 2 , 


2 _ ,2 „2 “e /I f v sin ( ' nT SN r ^ 


°SN " 2(KT s + VV B c'o ( r"7 5 f \2 

C B c (7rT SH f) 


df 


(13). 


For B T > 100, equation (13) reduces to 

0 bri 


4 b = (KT = + W* b = ri 


C SB 


= [(KT B ) 2 + 2KT B F_ + P 2 ] 

sc' s c R P J B T 


(lb) 

(15) 


c SB 


The mean noise power P^, is P^B^ (i.e., KTJB.J, and the mean signal plus 


D c 


noise power is KT s B c + P^. The mean signal power is P g ^ - P^ “ P R . 

Using the symbol N for the mean noise power, KT s B c , and the symbol S for 
the mean signal power, P^ equation (15) becomes 


<4 = <H 2 + 2HS + S 2 ) j-i 


(16) 


c SB 


In terms of the K/S ratio, equation (l6) can be written as 
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1 


(IT). 


Og K = S 2 [(K/S) 2 + 2 N/S + l] — 


c SH 


Similarly the equation for the variance of the noise cr^, is 


2 m 2 / 1 ' 
N c T K 


In terms of N/S this equation becomes 


4 - s2 t« s > 2 re 1 

c N 


(18). 


(19). 


Equations (l 6) and (l8) can be found in another way by defining the transfer 

function of the integrator by its equivalent bandwidth. As shown in figure 

A2-2,the equivalent bandwidth of the integrator is 1 /t, For B T > 100 the 

c 

power density P^B^ of the input signal to the integrator can be assumed 

constant across this equivalent bandwidth. The output power of the integrator 

2 2 

(i.e. the variance of the signal O ) is P B /t. Substituting in the signal 

u c 

plus noise power density for P , equation (16) is obtained, and substituting 
in the noise power density for P^, equation (l8) is obtained. 

The standard deviation of the measurement P^, can now be obtained by 
substituting equations (lj) and (19) into equation (7); thus, 

2 1/2 

G = S[=-“— (N 2 /S 2 + 2 N/S + 1) + (20) - 

M VSN Vn S 2 


Dividing through by S we have the normalizing factor, K , which gives the 
standard deviation of the estimate of P„ as a proportion of the true value 

A 

of P„. That is. 


K p S B t 


c SH 


(~ + 2 N/S + 1) 
S 


Vn S 2 


( 21 ). 
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so that the term 


The integration time T , will be some multiple of 
Kj. = T u/ T gu can ' be introduced in equation (21); i.e., 



1 


B c T SW 


(1 + 2 N/S + \ [1 + “-]) 
S 2 K t 


( 22 ). 


In general, as the integration time for noise is made longer, the 

standard deviation of the measurement will become smaller. The total integra- 
tion time for the signal plus noise is a function of the transmitted pulse 
period T , pulse repetition period T, and the measurement period t , and 

Jr P 

the noise measurement is performed during the off period between transmission 
pulses (see figure 27)- The signal plus noise integration period Tg^, will 
become smaller as is increased, thus increasing the standard deviation of 

the measurement. There is, however, an optimum choice of which will 

minimize the required S/N (maximum N/S) for any given measurement accuracy, 

K . By solving equation (22) for N/S, i.e.. 


N/s - 


1/2 

) 


K. 


K t + 1 


(23), 


and substituting in the equation for ig^ (from figure 27) 


SN 2R 


cmax 


t P T P 


(2U), 


+ K t (T P + 1} + \ 


the value of which gives the maximum N/S can be determined. The 

conventional means of determining the optimum value of — i.e., setting 
the derivative of equation (23) equal to zero and solving for K — results 
in a high order equation that is difficult to solve. A simpler approach, to 
determining an optimum value of K^, is through a graphical plot of the 
numerical evaluation of equation (23)- Using a K of .5, and the baseline 

lr 

design values for the various times in equation (2*0, and the bandwidths for 
doppler cells 1, 5, and 10, the graphs shown in figure A2-3 were determined. 
From these curves we see that the peah N/S occurs at a of approximately 

1.75- Also, the values of N/S are relatively constant for IC. values from 

A2-6 



It may be desirable to improve the measurement accuracy at the expense of 
resolution by averaging the measurements from multiple data cells. For 
example, this may be done In the high resolution design by averaging two adja- 
cent 25 tan resolution data cells to obtain a 50 km resolution measurement that 
has a standard deviation less than those for the 25 tan cells. This is shown 
as follows. Let 


and 



(25) 



( 26 ) 


be the normalized standard deviation of two adjacent cells. The errors 0.„ 

Ml 

and 0^ can be assumed to be independent. The average measurement of o° 
for the two adjacent cells is defined as, 


o o 
o °I +a 2 
°a - T~ 


(27). 


The standard deviation on this measurement will be 


0 ,. =1/2 (d(L + ofL) 

Ma Ml M2 


1/2 


(28): 


Substituting from equations (25) and (26), equation (28) becomes 


♦ ( Kp ^) 2 ] 1 ' 2 


(29). 


The normalized standard deviation, K , will be 



( 30 ) 


Ma 

XT — - 

HPa o 

[(^ g ^) 2 ♦ ( Kp ^) 2 ] 172 

«$ + 4 


(31). 


Assuming C 


c° (which is generally valid). 


K = 1/2 (K^ + \<? ) 
pa pi p2 


1/2 


(32), 


If K = K = K , then 
pl P2 P 


K 

K " “Z 

pa J2 


(33), 


For averaging a multiple number of data cell measurements, equation (32) can 
he written as 


K = 1/n (K" 
pa pi 


+ K + K _ 
p2 p3 


K2 ) 

pn 


1/2 


(3*0. 


If all the K 's in equation (3*0 are equal, which would generally he the case 
for measurements from contiguous data cells at a constant incidence angle, 
then equation (3*0 becomes 




/n 


( 35 ). 
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Ant. 

T a \7 


6 L R> T o 


L 


Loss 


T s 


TT 

T T 
out s 


F, T 


R 


Pre-amp. 

receiver 


T a = ant. noise temp. 

T q = amb. temp, of loss 

L„ = loss factor = ■ in P ut 
R output 

Tp = pre-amp-rec input noise temp. 

F = pre-amp-rec noise figure 

T = system noise temp, referred 
to pre-amp input 

T* = system noise temp, referred 
to ant. terminals 


Case I - Equation for T 


T s ■ T out + T R 


( 1 ) 


T ° Ut = X R (1 ' V T ° 
t r » T 0 < F - 1J 

+(1 -r- )T o + T o (F - ') (4 > 
R R 


(2) — Output noise temp, of a loss with input 
termination having noise temp., T - see 
ref. 32 *" a 




(5) - Loss Lp must be included in link 
calculation 


Case II - Equation for T* 


T’ = T + T 
s a e 


T e ■ < L R - D T o + L R T 


R 


T s - T A + < L R ‘ Wo + L R T R 


( 6 ) 


(7) - Input noise temp, to loss cascaded 

with amplifier - see ref. 32 

( 8 ) 


since T R = T Q (F - 1) 


Then T' = T + T (L„F - 1) 
s a o R ' 


— Loss Lp not included in link 
calculation 


Figure A2-1.- Derivation of system noise temperature equation. 




Input 


Sq. law 
detector 




Output 




Radar return signal and receiver noise assumed as Gaussian signals, 
uniformly distributed in I.F. BW B c 

For signal plus noise P D = KT g + P R /B c 

For noise only = KT g 

P R - P SN - P N * KT s B c + P R - KTB c “ P R 


Figure A2-2.- Spectrum at output of square law detector and gated integrator. 




Noise to signal ratio, N/S, db 
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Figure A2-3.- N/S requirement as a function of r^/rg^ ratio. 
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APPENDIX 3 

FAN BEAM ANTENNA 

By Dr. M. C. Bailey 
Langley Research Center 

This appendix presents the configuration and characteristics of a fan 
beam antenna, with application to the SeaSat A Satellite Scatterometer. 

DISCUSSION 


The basic antenna configuration is illustrated in figure A3-1. It 
consists of an array of pyramidal horns which are corrugated on two sides 
in order to produce equal beamwidths in the plane transverse to the ar‘ .y 
(ref. A3“I). For Mode I operation, the transverse dimension of the horn 
aperture (5.6 cm at 13.9 GHz) is chosen so as to produce a 25° half-power 
beamwidth. The array would have to be about 2.7 meters long (107 horns) at 
13.9 GHz in order to obtain a 0.5° half-power beamwidth in the longitudinal 
plane with -^30 dB sidelobes. 

The horn array elements are coupled to a square feed waveguide through 
centered crossed slots. When the feed waveguide is excited by orthogonally 
polarized waves (through a waveguide orthomode coupler) , the two arms of the 
cross slot provide the isolation and excitation for the radiating horn dual 
polarization. The tapered array illumination (for sidelobe reduction) is 
achieved by varying the slot size (slot-to-waveguide coupling (ref. A3-2)) 
along the feed waveguide. 

The computed radiation pattern in the transverse (yz) plane for the 
array of figure A3-1 is shown in figure A3-2 for the aperture electric field 
polarized along the x-axis (longitudinal to the array) , This pattern is that 
of an H-plane flared horn (ref. A3-3) since the corrugations only affect 
the electric field which is polarized normal to the corrugated surface. With 
proper design of the corrugations (ref. A3-4) , the ys-plane pattern for the 
other polarization should he very similar. 
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The array pattern in the longitudinal (xy) plane is determined primarily 

by the number of elements (N , the horn spacing (d) , the phase difference 

between adjacent elements (A), the element excitation distribution (C^) » and 

the element pattern (F($)) as 

H 1 
el 

E(§) = F(#) l C exp(j(n-l) (A-2iTd/A sin 0)) (1) 

- n 
n-1 


where $ is the angle measured off the z-axis and A is the free space 
wavelength. 

The radiation pattern is a maximum when the phase difference between 
the far-field signals of adjacent elements is an integral number of wave- 
lengths , i . e . , 


(A-2iTd/A sin $ ) = 2Trm m = 0,+l,+2, . . . (2) 

m — — 

where is the direction of the corresponding field maximum. 

In a slotted waveguide array, the phase difference between adjacent 
elements is determined by the slot spacing (s) and the guide wavelength 

(A ) , i . e . , 


A = 2Tis/A g . (3) 

Since the slot spacing is equal to the horn spacing (s = d) , substituting 
(3) into (2) and rearranging gives 

sin $ - A/A - mA/d. (4) 

m g 

There may be more than one value of m which will yield real values 
for § , corresponding to different beam maxima or "grating lobes. *' The 
position of the m = 0 lobe Is independent of the element spacing as seen 
from (4) with m = 0, i.e., 


sin § '« A/A . 
° g 


( 5 ) 
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However, within the normal bandwidth of rectangular or square waveguide 
(0.5 <_ A/A 0.866), the beam direction for the m = 0 lobe is restricted to 
the range of 30 to 60 from broadside. A broadside beam can be obtained, 
however, by selecting a slot spacing such that sin ^ = 0 for one of the 
other lobes. Then (4) gives, for a broadside beam, 

d = mA • (6) 

s 

Thus, by choosing a slot spacing equal to the guide wavelength of the feed 
waveguide, the peak of the m = 1 lobe will be in the desired direction. This 
is the criterion used for the design of the feed waveguide dimension and slot 
spacing of the array in figure 1. The value of A/A was chosen so as to 

°o 

separate the m = 0 lobe as far as practical ($ q —59 ) from the broadside 
beam. 

The xz-plane pattern for the array configuration of figure A3-1 with 
107 elements and with a -30 dB Taylor distribution (ref. A3-5) is shown in 
figure A3-3. This pattern was computed from (1) with A = 2ir and d = 1.165A. 
In addition to the m = 0 and m = 1 lobes, the m — 2 lobe is also present 
at an angle of 59° from broadside, toward the feed end of the array. The 
level of the lobes at 59° from broadside are determined by the horn element 
pattern at the corresponding angle. These lobes will be about -9 dB for the 
electric field polarized transverse to the array, and about -16 dB (or lower) 
for the orthogonal polarization. 

A rough estimate indicates that the return signal due to the "grating 
lobes" at 59° from broadside may be sufficiently reduced by the longer path 
length and smaller radar scattering cross-section of the surface so as not 
to cause any significant interference problem at an altitude of 808 km. 

An alternate configuration is illustrated in figure A3-4 in which the 
array of corrugated horns of figure A3-1 have been replaced by a corrugated 
trough. The radiation patterns should be very similar to those of the array 
of corrugated horns; however, the cost of fabricating a 2 to 3 meter long 
antenna may be appreciably reduced by selecting the corrugated trough design. 

For Mode II operation, which requires a 1/2° x 40° beamwidth, the basic 
antenna configuration can be modified by cutting off the horn elements (or 
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the trough for the alternate configuration) such, that the transverse aperture 
dimension is 3.95 centimeters instead of 5.6 centimeters. The tapered length 
then "becomes 5.2 centimeters instead of 11.0 centimeters. These dimensions 
are illustrated on figure A3-4 for the corrugated trough antenna. 

j 

i 
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Figure A3- 3- xz -plane radiation pattern of 107 element pyramidal horn array 
with -30 dB Taylor distribution and 1.1657^ spacing, 
(polarization in y-direction) 






